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I/"*) Abstract 
On 

J>~» The transition from reversible microdynamics to irreversible transport can be 

studied very efficiently with the help of the so-called projection method. We give a 
concise introduction to that method, illustrate its power by using it to analyze the 
00 well-known rate and quantum Boltzmann equations, and present, as a new applica- 

tion, the derivation of a source term which accounts for the spontaneous creation of 
ON electron-positron pairs in strong fields. Thereby we emphasize the fundamental im- 

portance of time scales: only if the various time scales exhibited by the dynamics are 
I/-} widely disparate, can the evolution of the slower degrees of freedom be described by 

a conventional Markovian transport equation; otherwise, one must account for finite 
memory effects. We show how the projection method can be employed to determine 
these time scales, and how -if necessary- it allows one to include memory effects 
-I— i in a straightforward manner. Finally, there is an appendix in which we discuss the 

i-^ concepts of entropy and macroscopic irreversibility. 
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1 Introduction 



Quantum transport equations confront us with both technical and 
conceptual problems. The macroscopic evolution of complicated quantum 
systems can often be described by means of transport equations. Prominent 
examples are the Langevin equation, Master equation, hydrodynamic equa- 
tions, rate equations, or the quantum Boltzmann equation. In this context the 
following technical questions arise: 

- Under which conditions are such transport equations applicable? Can one 
formulate precise mathematical criteria for their validity? 

- Can quantum transport equations be derived -in a rigorous fashion- from 
the underlying microscopic dynamics? And if so, how? 

Furthermore, the macroscopic evolution usually displays a high degree of ir- 
reversibility, which is reflected in a continuous increase of entropy. Associated 
with these features are the following conceptual problems: 

- What is the proper definition of the entropy? Why does its increase corre- 
spond to a physical thermalization of the system? 

- How can the irreversible behavior be reconciled with the reversibility of the 
underlying microscopic dynamics? 



In principle, these problems have long been solved. For the deriva- 
tion of quantum transport equations there are a variety of tools available, such 
as the projection method [1-17], the -closely related- Zubarev method [18- 
20], thermofield theory [21-23], Wigner function [24,25], and Green's function 
techniques [26-29]. Independent of the particular tool employed, one gener- 
ally arrives at the conventional transport equations only after certain approx- 
imations: coarse-graining, smoothing of rapidly oscillating functions, neglect- 
ing memory effects, to name just a few. All these approximations essentially 
amount to exploiting a separation of scales in the system [30,31]. Consequently, 
mathematical criteria for the validity of transport equations can usually be for- 
mulated in terms of these different scales. Such exploitation of widely disparate 
scales is not unique to quantum transport theory; it can already be found in 
classical transport theory. A famous example is furnished by the derivation of 
the classical Boltzmann equation: Boltzmann's "Stofizahlenansatz" (assump- 
tion that collisions may be considered statistically independent) crucially de- 
pends on the fact that the duration of each individual scattering process is 
much shorter than the average time that elapses between two successive col- 
lisions. There are hence two time scales, t co \\ and Tf ree , which must be widely 
disparate: t co \\ <C Tf ree . Not surprisingly, the simple Boltzmann description 
breaks down as soon as this separation condition is no longer satisfied, as is 
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the case in very dense or very strongly interacting plasmas. 

The separation of scales, in conjunction with a large number of degrees of 
freedom, is also responsible for the irreversible features of the macroscopic 
dynamics. The conceptual issues surrounding this macroscopic irreversibility 
are by now well understood; all relevant ideas have at some point been spelt 
out and can be found in various places in the literature [32-36]. 



Nevertheless, some of the issues have recently received renewed at- 
tention. It is impossible to present here an exhaustive overview over all 
the current activities in the field of quantum transport; yet our own area of 
specialization alone -the theoretical study of highly energetic nucleus-nucleus 
collisions- has recently witnessed many investigations of questions like the 
following: 

- Is the Boltzmann equation (or any other Boltzmann-like transport equation) 
applicable to high energy nucleus-nucleus collisions [37]? 

- More specifically, are there non-negligible memory effects [38-40]? If so, how 
do they affect the evolution of the system? Can they be accounted for in a 
non-Markovian generalization of the Boltzmann equation? How can such a 
non-Markovian Boltzmann equation be derived from first principles? 

- How can genuine quantum phenomena without classical analogues, such as 
spontaneous particle creation in a strong external field, be incorporated into 
existing kinetic equations [41-45]? How do they contribute to thermalization 
[46,47]? 

Aside from these technical questions, controversies continue over the origin and 
meaning of irreversibility [48]. Some of the recent debate has been stimulated 
by work on "decoherence" [49-56], i.e., the phenomenon that macroscopic 
systems obey essentially classical equations of motion despite the quantum 
mechanical nature of the underlying microscopic dynamics. 



One tool that appears especially suited for the study of such issues 
is the projection method. In order to be suitable for the investigation of 
the above issues, a theoretical framework should 

- permit an analysis of time scales; 

- if the scales are widely disparate: provide a well-defined, controlled pro- 
cedure to exploit this separation and to recover conventional Markovian 
transport theories; 

- if the scales are not widely disparate: allow for the derivation of non- 
Markovian transport equations and hence for the inclusion of finite memory 
effects; 
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- furnish a well-defined procedure to systematically expand any such non- 
Mar kovian correction in powers of time scale ratios; and 

- permit the definition of a relevant entropy and thus the quantitative study 
of irreversible behavior. 

One framework that satisfies all these conditions is the so-called projection 
method. As its name indicates, it is based on projecting the evolution of a 
quantum system onto some lower-dimensional subspace of the space of ob- 
servables (Liouville space). The projection method is broadly applicable and 
affords a systematic way of deriving quantum transport equations from first 
principles. With its help, one may obtain many seemingly unrelated equations 
such as the Langevin-, Master-, Boltzmann-, or time-dependent Hartree-Fock 
equation. The approach permits a systematic analysis of time scales and thus 
furnishes the desired criteria for the validity of the Mar kovian limit. Fur- 
thermore, the projection method allows for the incorporation of finite mem- 
ory effects and thus constitutes a useful tool for the study of non-Markovian 
corrections. Since it is formulated in the language of density matrices, mathe- 
matical manipulations and approximations within this formalism have a direct 
probabilistic, and therefore physical, interpretation, providing useful guidance 
during complicated calculations. Crucial concepts like the level of description, 
memory time, relevant entropy, etc., can be easily defined, making the projec- 
tion method especially suited for the issues surrounding the irreversibility of 
macroscopic processes. 



Our review will focus on specific aspects of the projection method. 

The projection method was originally developed by Nakajima [1], Zwanzig 
[2-6], Mori [7,8] and Robertson [9-13]. Its basic concepts have been reviewed 
[14] and are also presented in some textbooks on statistical mechanics [15-17]. 
Despite this wealth of available material we think that there are important rea- 
sons why another review might be called for: While in some branches of physics 
such as solid state physics [57] or quantum optics [58] the projection method 
has been employed extensively, there are other areas such as the physics of 
high energy nuclear collisions [59] in which quantum transport theory is used 
to describe the nonequilibrium dynamics of many-body systems, yet in which 
the projection method has found little or no application. In parts of the com- 
munity, therefore, the projection method and its potential advantages have 
remained largely unknown. 

For the reader who is not familiar with the projection technique we intend to 
give a concise introduction which focuses on the basic concepts, yet contains 
all essential formulae needed for applications. In order to demonstrate the 
power of the method we will present derivations -from first principles- of two 
well-known quantum transport equations, the rate equation for elastic scatter- 
ing and the quantum Boltzmann equation. We will further illustrate the use 
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of the projection method in linear response theory by deriving Drude's for- 
mula for the electrical conductivity. Aside from this general introduction and 
illustration we shall place particular emphasis on the following three aspects: 

(i) Time scales. We will discuss 

- why time scales play a fundamental role in the description of macro- 
scopic dynamics; 

- how the projection method can be used to determine these time scales 
(time scale analysis); 

- how, in this fashion, one can obtain quantitative criteria for the validity 
of conventional Markovian transport theories; 

- how -in cases in which these criteria are not fulfilled- the projection 
method allows one to derive generalized transport equations which are 
no longer Markovian and hence account for finite memory effects; and 
finally, 

- how such non-Markovian corrections may be expanded systematically 
in powers of time scale ratios. 

These aspects will be illustrated in the cases of the rate and quantum 
Boltzmann equations. 

(ii) New application. We will derive a (generally non-Markovian) source term 
in the quantum Boltzmann equation which accounts for the spontaneous 
creation of particle-antiparticle pairs in strong fields. 

(iii) Entropy and irreversibility. Even though -as we have claimed- all relevant 
ideas are already on the market, they are rarely presented in a coherent 
fashion. We wish to offer our own view of the subject, and discuss how 
the generally problematic issues (definition of relevant entropy, origin of 
irreversibility, etc.) present themselves in the context of the projection 
approach. In particular, we will emphasize the fundamental importance 
of disparate time scales and elaborate on the close connection between 
irreversibility and reproducibility of macroscopic experiments. 

Our goal is to elucidate these specific areas of interest, rather than to attempt 
a complete account of all actual or potential applications of the projection 
method (a task which would appear quite Herculean). For this reason we offer 
our apologies in advance to those numerous practitioners of the projection 
method who might not find their work adequately represented. 



The review is divided into five parts: 

(i) Following this introduction we first collect some theoretical preliminaries, 
notably the concepts of Liouville space, superoperators and generalized 
canonical states (section 2). 

(ii) Section 3 is then devoted to the essential ideas and general formalism 
of the projection method. We introduce the concept of projection su- 
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peroperators and derive the general projected equation of motion for an 
arbitrary set of selected observables. We discuss two particularly impor- 
tant realizations of this projection which give rise to, respectively, the 
Langevin-Mori and Robertson equations. We also include a discussion of 
possible approximations, with special emphasis on the Markovian limit. 

(iii) Next, in section 4, we turn our attention to specific examples and appli- 
cations. After some general remarks on quantum transport equations we 
present detailed investigations of the acceleration term in the quantum 
Boltzmann equation, the rate equation for elastic scattering, the Drude 
formula for the electrical conductivity, the source term in the quantum 
Boltzmann equation which accounts for spontaneous particle creation, 
and the collision term in the quantum Boltzmann equation. In all cases 
we proceed as follows: we first employ the projection method to derive the 
exact non-Markovian equation of motion; then perform a thorough time 
scale analysis; obtain criteria for the validity of the Markovian limit; and 
only then, where appropriate, recover the conventional Markovian trans- 
port equation. 

(iv) In section 5 we briefly summarize our findings and present an outlook on 
further potential applications of the projection method. 

(v) There is an appendix devoted to conceptual issues. In that appendix 
we begin with a general investigation of the entropy concept, and discuss 
why the entropy represents an appropriate measure to describe a system's 
approach to equilibrium. We shall argue that entropy -along with the so- 
called "relative entropy"- serves mainly as a diagnostic tool, a claim that 
we support with an account of various uses of the entropy concept in 
statistical physics. Only then we proceed to discuss temporal variations 
of the entropy, in particular the second law of thermodynamics and the H- 
theorem. We emphasize the intimate connection between irreversibility, 
reproducibility, and the existence of well-separated time scales. 
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2 Preliminaries 



2.1 Observables and Their Dynamics 



Observables are vectors in Liouville space. Linear operators acting on 
Hilbert space can be added and multiplied by complex scalars; they may thus 
be regarded as elements of an abstract vector space. This vector space of linear 
operators is called Liouville space [17]. Its Hermitian elements, i.e., Hermitian 
operators A, B } . . ., constitute observables. The state at time t of a quantum 
system is described by the statistical operator, or density matrix, pit) with 
properties 

P ] = P , P>0 , trp = l . (1) 

In the state pit) an observable A has the expectation value 

(A)(t) = tr[p(t)A] . (2) 



Liouville space is endowed with two scalar products: the Liouville space inner 
product, 

(A\B) :=tr [A^B] , (3) 



and the canonical (or symmetrized) correlation function [60] with respect to 
the state /?, 

l 

(A;B) P := J dXtx [p^tfp^B] . (4) 



The canonical correlation function with respect to a canonical equilibrium 
state p = Z~ x exp(—(3H) is also known as the Mori product. In terms of these 
scalar products, expectation values can be written as 

(A) = (p\A) = (l;A) p . (5) 



The dynamics is determined by the Liouville-von Neumann equa- 
tion. The dynamics of a quantum system is governed by its Hamiltonian 
H. In the Schrodinger picture the statistical operator evolves according to the 
Liouville-von Neumann equation 

^f = 4M • (6) 
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This equation has the formal solution 
p(t) = U(t,t )p(t )U\t,t ) 



(7) 



with U (t, t ) being the unitary evolution operator. A state p is called stationary 
if [H } p] = 0; while an observable A is called a constant of the motion if 
[H, A] = 0. Unless explicitly stated otherwise, we will always assume that 
neither the Hamiltonian nor the observables are explicitly time-dependent. 



The Liouville-von Neumann equation can be formulated and solved 
in the language of superoperators. An operator O acting linearly on 
vectors in Liouville space, A — > OA } is called a super operator. Three important 
examples are the Liouvillian C, the evolution superoperator U and the Green's 
function superoperator Q. The Liouvillian is defined as 

CA:=h- l [H,A] . (8) 



With its help the Liouville-von Neumann equation can be written in the form 



while the equation of motion for expectation values takes the form 

±(A) = i( P \CA) . (10) 

The Liouvillian is Hermitian with respect to the Liouville space inner product 
(•|-), but generally not with respect to the canonical correlation function (•; ■) p 
(unless p is stationary). 

The evolution superoperator is defined as 

U(t,t )p(t ) := U(t,t )p(t )U\t,t ) . (11) 

It determines -at least formally- the evolution of expectation values via 

(A)(t) = (p(t )\U(t ,t)A) . (12) 

The evolution superoperator is unitary with respect to the Liouville space 
inner product, but in general not with respect to the canonical correlation 
function. It has the property 

U- 1 (t,t )=U(t ,t) , (13) 
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and it satisfies the equation of motion 

j^(Mo) = -iJCU(t,t ) (14) 



or 



^U(t ,t) = iU(t ,t)C , (15) 



respectively, with initial condition U(t 0} t ) = 1. Multiplication with a step 
function yields the causal evolution superoperator 

W<(* ,<):=W(*o,<M(*-*o) (16) 



(where '<' means H < t r ). 

Finally, the (causal) Green's function superoperator is defined as the Fourier 
transform of the ( 1) evolution superoperator, 



^ (< )(t ,cu) := J dt exp(zcut)W (<) (t ,t) . (17) 

— oo 

Provided the Liouvillian C is time-independent and hence 

U{t ,t) = exp[i£-(t-t )] , (18) 

the associated Green's operators are given by 

£ (< )(t ,cu) = exp(iujto) G(<)(u) (19) 



wi 



ith 



g(u):=2iTS(u + C) 
i 



w + if] + C 



(20) 



>o+ 



Causal and ordinary Green's function superoperators are related via 



g < ( u ) = ±g( u ) + iP I (-±-^) , (21) 
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where 'Pr' denotes the principal value. With respect to the Liouville space 
inner product or any equilibrium correlation function, is Hermitian, while 
[i Pr . . .] is anti-Hermitian. 

2.2 Generalized Canonical States 

Generalized canonical states contain information about selected ob- 
servables only. At any given time t let the expectation values of some 
selected observables {G a } } 

g a (t):=(G a )(t) , (22) 

be the only information available. This information is generally not sufficient 
to reconstruct unambiguously the statistical operator p(t); there are many pos- 
sible choices. However, one choice of particular importance is the generalized 
canonical statistical operator 

Pcw [9a(t)]:= Z(t)- l e W {-^(t)Ga) ■ (23) 

Here Zit) denotes the partition function 

Z(t) := ti exp(-\ a (t)G a ) , (24) 

a summation over a is implied (Einstein convention), and the time-dependent 
Lagrange parameters {A a (t)} are adjusted so as to yield, at each time t, the 
correct expectation values for the selected observables. 

Among all states which satisfy the constraints tip = 1 and ti(pG a ) = g a} 
the generalized canonical state is the one which maximizes the von Neumann 
entropy 

S[p] ■= -ktr(plnp) . (25) 

For this reason the generalized canonical state may be considered "least bi- 
ased" or "maximally non-committal" with regard to missing information. A 
more detailed discussion of this maximum entropy rationale can be found in 
the appendix. 

Expectation values, Lagrange parameters and entropy satisfy ther- 
modynamic identities. The following relations hold at each time t. The 



10 



expectation values {g a } of the selected observables can be obtained from the 
partition function via 

g^-^-JnZ . (26) 



Inserting the generalized canonical state into the definition of the entropy 
yields 

S = klnZ + k\ a g a , (27) 



with differential 

dS = k\ a dg a . (28) 

Arbitrary expectation values (A) in the generalized canonical state /? can depend 
on the Lagrange parameters {A a }. If these vary infinitesimally, the expectation 
value (A) changes according to 

d(A) = -(SG a ;A)d\ a . (29) 

Here the correlation function (•; •) is evaluated in the state p ca ,m and the ob- 
servable 8G a is defined as 

SG a := G a - (G a ) . (30) 

With the help of the correlation matrix 

C ab := (8G a ;8G b ) = --^ (31) 

which relates infinitesimal variations of A and g } 

dg b = -d\ a C ab , d\ a = -dg b (C- l ) ha , (32) 

one may expand the expectation value also in dg: 

d(A) = (C- 1 ) ba (5G a ;A)dg b . (33) 
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3 Projection Method 



3.1 Motivation 

The dynamics is projected onto the level of description. When deal- 
ing with a complicated quantum system such as a gas or plasma, one is typi- 
cally neither interested in, nor capable of, describing the time evolution of all 
its microscopic properties. Rather, two typical problems are 

(i) in linear response theory: to determine dynamical susceptibilities, which, 
according to Kubo's formula [61,62], amounts to calculating dynamical 
correlations of the form {G a ;Q < (uj)Gb}eqi and 

(ii) more generally, away from equilibrium: to determine the evolution of 
certain selected expectation values {(G a )(t)}. 

Both calculations require knowledge of the dynamics of only a small set of 
selected ("relevant") observables {G a }. These relevant observables, together 
with the unit operator, span a (relatively small) subspace of Liouville space, 
the so-called level of description. 

The full statistical operator generally contains a lot of information which is 
not at all related to the relevant observables. This irrelevant information com- 
plicates calculations without providing any insight into the dynamics of the 
relevant observables. To discard such unnecessary "baggage," one projects the 
motion of the system onto the level of description. This projection yields a 
closed equation of motion for the relevant observables only; which, however, 
is generally no longer Markovian. 



The projection method allows one to exploit a separation of time 
scales. Typically, neither the original Liouville- von Neumann equation nor 
the projected, non-Markovian equation of motion for the relevant observables 
can be solved exactly. One must resort to clever approximations, a task which, 
in practice, amounts to finding suitable expansion parameters. Two possibili- 
ties come to mind: 

(i) The Liouvillian may have the form C = C° + aV with V denoting a per- 
turbation and a a small coupling constant. In this case one may expand 
in powers of a (perturbation expansion). 

(ii) Macroscopic systems may exhibit widely separated time scales r s h or t and 
Tiong- In this case an additional expansion parameter is furnished by the 
ratio (Tshort/'Hong)- As will be shown later, the zero-th order of such an 
expansion corresponds to the Markovian limit. 
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Table 1 

Examples for sets of relevant observables, associated relevant parts of the statistical 



operator, and 


rinds of discarded information 


Physical 
system 


Relevant 

observables {G a } 


Relevant part p re \(t) 


Discarded 
information 


^x&x. . . 

several 
subsystems 


all S 8 -observables 
(i = 1,2,...) 


p 1 (t)®p 2 (t)®... 


correlations 

between 

subsystems 


many- 
particle 
system 


all or selected 
single-particle 
observables 


r / 1 / ^ \ \ ? <} / 1 \ T \ 

Z exp(-X)A J - 7 (i)ajaj) 
(Hartree form) 


particle 
correlations 


S xB 
system in 
heat bath 


1 1 Ci 1 11 

all b -observables 
and bath energy 


9 / _l\ ^ r / — 1 / n / j_\ TTh\ 

p (t)(g)Z 1 exp(-p(t)H ) 


j_ 11 

internal degrees 
of freedom of the 
heat bath; 
system-bath 
correlations 


arbitrary 


all observables 
commuting with 
A = £i a t P t 




coherence with 
respect to A 



While the original Liouville-von Neumann equation may lend itself to a per- 
turbative treatment, it is not a good starting point for an expansion in powers 
of time scale ratios. In contrast, the projected equation of motion -with the 
influence of irrelevant degrees of freedom mapped onto a non-local behavior 
in time- allows one to detect and systematically exploit a separation of time 
scales. This is the basic practical merit of the projection method. 



3.2 General Formalism 

3.2.1 Relevant Part of the Statistical Operator 

Discarding irrelevant information yields the relevant part of the sta- 
tistical operator. The relevant observables have time-dependent expecta- 
tion values g a (t) = (G a )(t). Given these expectation values as constraints, one 
may at all times t define an associated canonical state p ca ,n[ga(t)]. Due to the 
entropy maximization involved, it contains information pertaining to the rele- 
vant observables only; irrelevant information has been discarded entirely. For 
this reason the generalized canonical state is also called the relevant part of 
the statistical operator and denoted by p re i(t). Some common examples for sets 
of relevant observables, associated relevant parts of the statistical operator, 
and kinds of discarded information, are listed in table 1. 
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Fig. 1. Arbitrary observables are projected onto the level of description. The pro- 
jection may vary in time. 

3.2.2 Projectors 

Arbitrary observables are projected onto the level of description. 

A projector in Liouville space is a superoperator V satisfying V 2 = V. (It 
need not be Hermitian with respect to any scalar product.) If V is a projector 
then so is its complement Q := 1 — V. For our purposes we consider projectors 
which project arbitrary vectors in Liouville space onto the level of description; 
i.e., for which 



The relevant observables and therefore the level of description are assumed 
to have no explicit time-dependence. In contrast, one allows the projector to 
depend on the expectation values {g a (t)} of the relevant observables, thus 
making it an implicit function of time: V = V[g a (t)] (cf. Fig. 1). Since the 
level of description itself remains fixed, a second projection will have no effect, 
even at a different time: 



Two important examples are the Kawasaki-Gunton and the Mori 
projector. In principle, there are many possible choices for the projector, 
but the following two are particularly useful: 

- The Kawasaki-Gunton projector [63] is defined such that the projection 
onto the level of description is orthogonal with respect to the canonical cor- 
relation function (•; -) re i(t). The correlation function is defined in the time- 
dependent relevant part of the state, p re i(t). The Kawasaki-Gunton projector 



VA = A iff A G span{l, G a } 



(34) 



V[g a {t2)]V[g a {t l )] = V[g a {t l )} V t u t 2 



(35) 
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thus varies in time: it depends -via p Te \(t)- on the time-dependent expecta- 
tion values {g a (t)}. Explicitly it reads 

V KG [9a(t)] A := (1; A) reiw l + (C^iSC^^AUuSOF® . (36) 

Here the operators SG™^ are defined as 

8G^) :=G a -g a (t) , (37) 

and C re i(t) denotes the correlation matrix evaluated in p re i(t). 
- The Mori projector [7] is defined such that the projection onto the level of 
description is orthogonal with respect to the canonical correlation function 
(•; -) eq . The correlation function is now defined in equilibrium, p eq . The Mori 
projector is thus time-independent. Explicitly it reads 

V M A := (1; A) eq l + {C- l ) ba {8G a] A) eq SG b . (38) 
Here the operators 8G a are defined as 

SG a := G a - (G a ) e(i , (39) 

and C denotes the correlation matrix evaluated in equilibrium. 

The above two projectors ('KG' and 'M') lead to different expectation values 
of the projected observable: 

f (^)reim : KG 

(p(t)\VA)=\ K /relW (40) 

( (A) eq + (C- v ) ba (SG a ; A) eq ■ [g b (t) - (G b ) eq ] : M 

Provided the equilibrium state is completely characterized by the expectation 
values {{G a ) eq } } i. e., provided the relevant observables contain all constants 
of the motion as a subset, then according to (33) the two expectation values 
agree to first order in [gb(t) — {Gb) e q] - In other words, in the linear regime close 
to equilibrium the Kawasaki-Gunton and Mori projectors yield identical re- 
sults. Far from equilibrium, however, the two projections may have completely 
different effects. Which of the two projectors is better suited to a particular 
application differs from case to case; their respective advantages and disad- 
vantages will be discussed in section 3.4.2. 

3.2.3 Equation of Motion for the Relevant Observables 

Eliminating irrelevant degrees of freedom gives rise to a memory 
term and a residual force term. We return to the generic case of an 
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arbitrary, possibly time-dependent projector V(t) = V[g a (t)]. The only specific 
property that we require is that it satisfy 

(p(t)\±V(t)A) = VA . (41) 



Both the Kawasaki-Gunton and the Mori projector fulfill this requirement. 
As before, we assume that the Hamiltonian and hence the Liouvillian are not 
explicitly time-dependent. 

While U(t ,t) is the evolution superoperator corresponding to the Liouvillian 
£, let T(t ,t) be the evolution superoperator corresponding to 

Arr(<) := Q(t)£Q(t) ; (42) 



i. e., it satisfies the equation of motion 
d 

—T{t ,t) = iT{to,t)£ irr {t) (43) 



or 



^- i T(t\t) = -iCUt')T(t\t) , (44) 



respectively, with initial condition T(t ,t ) = 1. Since Q projects out the 
irrelevant component of an observable, the superoperators C irr and hence T 
may be pictured as describing the evolution of the irrelevant degrees of freedom 
of the system. Using the definition (44) of T, the Liouville-von Neumann 
equation (9), as well as the property (41) of projectors, we find 



( p (t)\Q{t) = {p{t )\Q{t )T{t< h t) 

^(p(t')\Q(f)) T(t',t) + ( P (t')\Q(t')-^T(t',t) 



+ I dt ' 



to 

t 

= (p(t )\Q(t )T(t ,t) + i J dt'(p(t')\V(t')jCQ(t')T(t',t) . (45) 

to 

Without changing the left-hand side we may multiply this equation by Q(t) 
from the right. The equation of motion for the expectation values of the rele- 
vant observables can then be written as 



m={ t (p(t)\G«) 
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= i(p(t)\V(t)CG a ) + i(p(t)\Q(t)CG a ) 

t 

= i( P (t)\V(t)CG a )- J ' M{ P {t')\V{t')CQ{t')T{t',t)Q{t)CG a ) 

to 

+ i(p(t )\Q(t )T(t ,t)Q(t)CG a ) . (46) 

This is the desired result. The equation of motion is exact: no physical as- 
sumptions or approximations have been made. We have largely succeeded in 
replacing (p\ — > (p\V, thus eliminating almost completely the influence of irrel- 
evant degrees of freedom. The only exception is the last term which accounts 
for the residual influence of irrelevant components in the initial state; but in 
many applications this term can be shown to vanish, or else to be manageable. 

Of course, the elimination of irrelevant degrees of freedom has its price. Com- 
paring the above equation of motion with the Liouville-von Neumann equa- 
tion, we notice that it contains two terms which are qualitatively new: (i) an 
integral ("memory") term, containing contributions from all times between 
the initial and the present time; and (ii) a "residual force" term describing 
the effect of irrelevant components in the initial state. The physical meaning 
of these two terms can be easily discerned if read from left to right: (i) At 
time t' < t relevant degrees of freedom (projected out by V) couple via an 
interaction (C) to irrelevant degrees of freedom (projected out by Q), which 
subsequently evolve in time (T) and, due to a second interaction (£), acquire 
relevancy again, thus influencing the evolution of the relevant observable G a 
at the present time t. (ii) Irrelevant components in the initial state (Q) evolve 
in time (T) and, due to interaction (£), acquire relevancy at the present time 
t. 



The projected equation of motion is a good starting point for ap- 
proximations. In essence, the influence of irrelevant degrees of freedom has 
been mapped onto a non-locality in time: future expectation values of the rel- 
evant observables are predicted not just on the basis of their present values, 
but based on their entire history. At first sight such an equation of motion 
may seem even more complicated and awkward than the original Liouville- 
von Neumann equation; and indeed, in many cases it is just as hard to solve 
exactly. However, mapping the complexity of the problem onto a non-local 
behavior in time opens the way to the systematic exploitation of a possible 
separation of time scales. Therefore, the above equation of motion serves as 
a much better starting point for various approximations, in particular the 
Markovian approximation. In fact, many of the well-known equations of non- 
equilibrium statistical mechanics can be derived from there, either as special 
cases or with the help of a Markovian approximation. Examples that we will 
discuss in detail are the Robertson and Langevin-Mori equations and, later 
on, the rate equation for elastic scattering and the quantum Boltzmann equa- 
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tion; other examples being the Nakajima-Zwanzig, Master, or time-dependent 
Hartree-Fock equations [17]. 

3.2.4 Approximations 

Two important approximations are perturbation theory and the 
Markovian limit. The projection method was initially motivated by the 
claim that it permits a new kind of approximation based on the existence of 
disparate time scales in the system: the Markovian approximation. The gen- 
eral equation of motion for the relevant observables thus serves as the starting 
point for two important approximations: 

- the conventional perturbation expansion; and 

- the Markovian limit, as well as (associated with it) the quasistationary limit. 

These two approximations will be discussed separately. In order to streamline 
notation, we will denote the kernel of the memory term by 

1C{t\t):=V{t')CQ{t')T{t\t)Q{t)CV{t) . (47) 



Perturbation theory is based on a suitable decomposition of the 

Liouvillian. Often the Liouvillian can be split into a free part and an in- 
teraction part, 

C = £(°) + V , (48) 

corresponding to a decomposition H = + V of the Hamiltonian. Provided 
free evolution does not mix relevant and irrelevant degrees of freedom, i.e., 
provided 

[£ {o \V(t)] = , (49) 

then in the memory term VCQ = VVQ and QCV = QW. In that case the 
memory term is at least of order (9(V 2 ); it can be dropped in first order pertur- 
bation theory. One prominent application of first order perturbation theory is 
to a system of nonrelativistic fermions with effective 4-fermion interaction: if 
one chooses the level of description to be spanned by all single-particle observ- 
ables and employs the Kawasaki-Gunton projector, then dropping the memory 
term will yield the time-dependent Hartree-Fock equations [64,65]. In this ex- 
ample first order perturbation theory amounts to a mean field approximation. 
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In order to obtain higher order contributions one must include the memory 
term and expand 

Q(t')T(t',t)Q(t) = Q(t')T {0) (t',t)Q(t) + 0(V) , (50) 

where is a time-ordered exponential of C[^(t). Using Q(t 2 )Q(ti) = Q(t 2 ), 
all the Q's appearing in 

T (o) 

can be shuffled to the left and absorbed into Qit' ), 
allowing one to replace simply by U^°\ To second order we thus obtain 

g^(t) = l (p(t)\T(t)CG a ) 

t-t 

- J dr (p(t-T)\V(t-T)VQ(t-T)uW(0,T)Q(t)VG a ) 



+ residual force . (51) 

The Markovian limit amounts to ignoring the system's past history. 

We have found that predictions of future expectation values of the relevant 
observables generally depend in a complicated manner on both their present 
values and their past history. There are thus two distinct time scales: (i) the 
scale T re i -or several scales {t^j}- on which the expectation values {g a {t)} 
evolve; and (ii) the memory time r mem which characterizes the length of the 
time interval that contributes significantly to the memory term. Loosely speak- 
ing, the memory time determines how far back into the past one has to reach 
in order to make predictions for the further evolution of the relevant observ- 
ables. Often the memory kernel )C(t',t) is proportional to some distribution 
(e. g., Gaussian or Lorentzian) in it — t'), in which case the memory time can 
be identified with the width of that distribution. Symbolically: 

£(*',*) «0 if (t-f)>T mem . (52) 



If this memory time is small compared to the typical time scale on which the 
relevant observables evolve, r mem <C r re i, then memory effects can be neglected 
and predictions for the relevant observables be based solely on their present 
values. One may then assume that in the memory term g a (t') ~ g a (t), and 
hence replace 

n9a(f)]^V[g a (t)} 
(p(t')\V(t')^( P (t)\V(t) 

T(t' } t)^exp[iQ(t)£Q(t)-(t-t')] . (53) 

With these replacements, the equation of motion -ignoring the residual force 
term- simplifies to 
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gJ m \t) = i(p(t)\V(t)£G a ) 

t-t 

-(p(t)\V(t)CQ(t) J dr exp[iQ{t)£Q{t) ■ r]Q(t)CG a ) . (54) 
o 

The change of expectation values of the relevant observables now depends 
only on their present values, their past history being completely ignored. This 
is the Markovian limit. Note that nevertheless the kernel K, still contributes: 
ignoring the system's past history does not mean that the memory term is 
eliminated altogether. If on a given level of description it is justified to take 
this Markovian limit, the level is called a Markovian level of description. 

Closely related to the Markovian limit is the quasi-stationary limit. For times 
t with it — t ) ^> T mem the exact "location" of the initial time t becomes 
irrelevant; one may then let t — > — oo. This removal of the cut-off t restores 
temporal homogeneity and hence energy conservation. 



One may systematically expand around the Markovian limit. Ex- 
panding in the memory term g a (t') around g a (t) } 

ga(t , )=g a (t) + g a (t)-(t-t , ) + ... , (55) 



will yield corrections to the Markovian equation of motion. These non-Marko- 
vian corrections are of successively higher order in (r mem /r re i). To illustrate 
this, we consider a simple model: the evolution of a single expectation value 
g(t) governed by a non-Markovian equation of motion 

t-t 

g(t) = - ( dT 7 (T)g(t-r) . (56) 



For simplicity we assume the kernel 7(7") to have the form 

r 

7( r ) = exp(-r/r me m) (57) 



with a finite memory time r mem . As long as r mem <^ it is justified to take 

the Markovian limit; one may then replace git — r) — > git) and obtain the 
approximate equation of motion 

g^(t) = -Tg(t) . (58) 
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If one wants to go beyond this approximation one must expand 



oo i ite 

9(t-r)=g(t) + ^(-T) k — kg (t) 



(59) 



which yields the exact equation 



m = 9 (m) (t)-vY J ( 



CO 




dt k 



(60) 



k=i 



Non-Markovian corrections are thus associated with successively higher (A;-th) 
order derivatives of g(t). Provided the expectation value evolves on a typical 
time scale r re i, it is 



3.3 Which observables should be relevant? 

The proper choice of relevant observables is dictated by the time 
scales. In principle, the choice of the set of "relevant observables" is arbi- 
trary. However, some choices are more useful than others. At first sight, it 
seems obvious that always those observables should be considered relevant 
which are actually being monitored in the experiment; but this choice is not 
always the most suitable. As an example consider an experiment performed on 
a dilute classical gas, in the course of which local energy and particle densities 
are being monitored. A theoretical description of this experiment will most 
likely start from the Boltzmann equation — even though in the Boltzmann 
equation all single-particle observables are taken to be relevant, not just those 
which are actually being monitored. The reason for employing the Boltzmann 
equation, rather than an equation of motion for the local densities only, lies in 
the fact that the level of description spanned by all single-particle observables 
is Markovian, while the smaller experimental level is not. 

Such a situation is quite common. Typically, the experimentally monitored 
observables are contained in some larger set of slowly varying observables. 
(For an observable to be measurable in practice it is usually necessary that 




(61) 



and hence the A;-th non-Markovian correction is of order 




(62) 
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it vary slowly.) While the experimental level of description is not Markovian, 
the level associated with the larger set is. In this case it is best to extend 
the level of description from the experimental to this larger Markovian level. 
Only this allows one to make use of the powerful Markovian approximation, 
a significant advantage in any explicit calculation. The proper choice of the 
relevant observables is thus dictated by the time scales. 

3.4 Representations of the Equation of Motion 

Representations of the general equation of motion can be obtained by choos- 
ing specific projectors V(t). We will discuss two choices: the Kawasaki-Gunton 
projector, which leads to the Robertson equation [9-12,63]; and the Mori pro- 
jector, which leads to the Langevin-Mori equation [7,8]. In both cases we will 
obtain the respective representation by inserting Eq. (40) into Eq. (46). 

3.4.I Robertson and Langevin-Mori Equations 

Inserting the Kawasaki-Gunton projector into the general equation 
of motion yields the Robertson equation: 

t 

g a (t) = i(£G a ) rem - J dt' (iCQ(t')T(t', t)Q(t)CG a ) rem 

to 

+(T(t (h t)Q(t)uCG a ) irr{to) . (63) 

Here 

Pirr{to) ■= p(to) ~ Pre\{to) (64) 

denotes the irrelevant part of the initial state. With the help of the identity 
{iCA) rem = {iCG a ;A) rem X a (t) VA , (65) 

and using the fact that the projector Q(t') is Hermitian with respect to 
(■; ')rei(t')? one may also write 

t 

g a (t) = (iCG c - G a ) rem \ c (t) + J dt' {Q(t')iCG c - Tit', t)Q{t)iHG a )^ {V) \ c {t r ) 

to 

+(T(t ,t)Q(t)i£G a ) irr{to) . (66) 
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Inserting the Mori projector into the general equation of motion 
yields the Langevin-Mori equation: 

t-t 

g a (t) = n b Jg b (t) - [ dT 7 b a (T)Sg b (t-T) + (T(t ,t)f a ) p{to) . (67) 



Here 

^:=-(C- 1 ) 6c (i£a;G a ) eq (68) 



denotes the frequency matrix; 

l h a (T):=(C- l ) hc (f c -T(^T)f a U (69) 

the memory matrix; 

f a := QlCG a (70) 

the stochastic force; and finally 8g a (t) is defined as 

Sg a (t) := g a (t) - (G a ) eq . (71) 

The Langevin-Mori equation and the second version (66) of the Robertson 
equation have similar structures. The former can be obtained from the latter 
by replacing 

(•; - W) -> (•; -) eq , Q(t)iCG a -> f a , X c (t) -+ -(C-^Sg^t) . (72) 

Langevin-Mori theory leads to an equation of motion for the dy- 
namical correlations. The Langevin-Mori equation may also be read as an 
operator equation if one replaces the expectation values 8g a (t) by Heisenberg- 
picture operators 8G a (t). The resulting operator equation can be inserted into 
the definition of the normalized dynamical correlation matrix [17] 

~ a (t):=(C-T(5G c -8G a (t)U , (73) 

with S^(0) = 8 ^, to obtain its equation of motion 

t 

S(f ) = !)-(*)- |dr 7 (r)S(t-r) . (74) 
o 
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Here we have employed a compact matrix notation. Note that both the initial 
time t and the residual force term have dropped out. 



Multiplying and ^(t) by step functions yields the causal normalized dy- 
namical correlation matrix S<(t) and causal memory matrix 7<(t), respec- 
tively. Their Fourier transforms are given by 



As a consequence of the equation of motion (74), the Fourier transformed 
S < (cu) can be expressed directly in terms of the frequency and memory ma- 
trices: 



In the Markovian limit the memory matrix falls off rapidly in time, hence its 
Fourier transform is practically constant; one may then replace 7<(cu) with 



3.Jf.2 Which representation should one choose? 

Whether the Kawasaki-Gunton or the Mori projector is more ad- 
vantageous depends on the specific application. Since the Robert- 
son equation and the Langevin-Mori equation are both exact -although in 
practice they can only be solved with suitable approximations- the choice of 
projector is not a matter of principle but of convenience. Depending upon 
the specific physical process under consideration, either the Kawasaki-Gunton 
projector and hence the Robertson equation, or the Mori projector and hence 
the Langevin-Mori equation, may serve as better starting points for explicit 
calculations. 

The Kawasaki-Gunton projector has one major disadvantage: 

- The Robertson equation leads to a system of coupled integro-differential 
equations for the expectation values {g a (t)}. Since the Kawasaki-Gunton 
projector is not constant in time, but depends in a non-trivial way on the 
{g a (t)} } these equations are generally nonlinear. 

On the other hand, employing the Kawasaki-Gunton projector has three po- 
tential advantages: 

- Often the initial macrostate can be characterized completely by the expec- 
tation values {{G a )(t )} of the relevant observables. The initial state is then 




(75) 
(76) 




(77) 



7<(0). 
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determined by maximizing its von Neumann entropy, and hence has just the 
generalized canonical form. In that case the initial state has no irrelevant 
component, which in turn implies that there is no residual force. Eliminating 
the residual force term in this manner greatly simplifies the calculation. 

- If all relevant observables commute, i.e., if [G a ,Cr&] = 0, then also the first 
term of the equation of motion vanishes: 

( P rei\CG a ) k ti(exp(- \ b G b )[H } G a ])=ti([G a} exp(-\ b G b )]H) = 0.(78) 

This in turn implies that VCV = and that hence in the memory term 
of the Robertson equation one may omit all complementary projectors Q. 
However, these conclusions must be taken with a grain of salt: their deriva- 
tion depends on the cyclic permutability of operators under the trace, a 
property which holds only if both G a and H are bounded [66]. Indeed, 
in section 4.2 we will encounter an example -the acceleration term in the 
quantum Boltzmann equation- where the above line of reasoning does not 
apply. 

- Many transport theories have as their relevant observables some set of num- 
ber operators {A^} or, more generally, of single-particle observables {a}ctk}. 
Evaluation of the memory term then requires calculating expectation values 
(a\a 3 a\ai . . .} re i(t') of products of several field operators. Because the relevant 
part p Te \(t') has the so-called Hartree form [67] (cf. table 1), Wick's theorem 
can be applied and these expectation values be evaluated. 

The advantages and disadvantages of the Kawasaki-Gunton projector are in a 
sense complementary to those of the Mori projector. The main disadvantage 
of the latter is: 

- The residual force term generally does not vanish, nor can it be calculated 
exactly. Rather, in most applications it must be suitably modelled. 

This has to be weighed against its advantages: 

- Due to the time-independence of the Mori projector, the Langevin-Mori 
equation has a convenient linear structure. This linearity permits Fourier 
transformation and thus a treatment of the dynamics in the frequency rep- 
resentation. 

- Frequency and memory matrices are defined in equilibrium and may there- 
fore be easier to evaluate. 

- In the context of linear response theory, evaluation of the Kubo formula 
can often be reduced to calculating the dynamical correlation matrix £<(u;) 
of some selected observables. This dynamical correlation matrix, in turn, is 
most easily calculated in the Mori formalism (Eq. (77)). There the prob- 
lematic residual force term no longer appears. 

Summary. To put it in a nutshell, there is usually a trade-off: complicated 
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features of the dynamics may appear either in the form of nonlinearities 
(Kawasaki-Gunton) or a non- vanishing residual force (Mori). Which projector 
one should choose depends on which complication is easier to handle. As a 
rule of thumb, the Kawasaki-Gunton projector is usually better suited for the 
description of processes far from equilibrium, while the Mori projector is bet- 
ter adapted to describing processes in the linear regime close to equilibrium, 
and to evaluating the Kubo formula. 

3.5 Accompanying Entropy 

Responsible for entropy generation is the memory term. Associated 
with the relevant part of the statistical operator is the relevant entropy 



This relevant entropy is also known as the accompanying entropy with respect 
to the level of description [17]. Its value generally varies in time. According to 
Eq. (28), relevant entropy is produced at the rate 



This formula can be evaluated explicitly by inserting the Robertson equation 
(66). We assume that the initial state is completely specified by the initial 
expectation values of the relevant observables so that there is no residual force. 
One can show that the first term in the Robertson equation never contributes 
to entropy generation; hence the only non-vanishing contribution stems from 
the memory term: 



This expression is valid arbitrarily far from equilibrium. 



In the Markovian limit the accompanying entropy increases mono- 
tonically. This statement is true generally; however, for simplicity, we will 
prove it here in the linear regime close to equilibrium. In that regime the above 
expression for the entropy production rate can be further simplified: one may 
replace Kawasaki-Gunton by Mori projectors and evaluate the canonical cor- 
relation function in equilibrium. Taking the quasistationary limit (t — > — oo), 



Srel[g a (t)] ■= -ktr(p re l(t)llip re l(t)) 



(79) 



Sr e i[g a (t)] = kX a (t)g a (t) 
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the entropy production rate then reads 

CO 

S^[ 9 a(t)] = k J drA c (t-r)(/ c ;r<(0,r)/ a ) eq A a (t) . (82) 

— CO 

Provided the evolution is Markovian, this reduces to 

S^\9a(t)] = 7rk{X c (t)f c ; 8(QCQ) X a (t)f a ) eq > . (83) 

The superoperator S(QCQ) and therefore the entropy production rate are 
non-negative, Q.E.D. 

This result is just one particular manifestation of the far more general H- 
theorem. The _ff-theorem, the second law of thermodynamics, and their con- 
ceptual ramifications will be discussed in greater detail in the appendix. 
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4 Application: Quantum Transport Equations 

4.1 Introduction 

When applying the projection method one generally proceeds in 
eight steps: 

(i) Choose the level of description by identifying the relevant observables. 

(ii) Determine, or make assumptions about, the initial state. 

(iii) Determine the microscopic dynamics (Hamiltonian). 

(iv) Choose the projector and hence the representation of the equation of 
motion. 

(v) Inserting the specific choice of relevant observables, initial state and 
Hamiltonian, obtain the exact, generally non-Markovian equation of mo- 
tion. 

(vi) First part of the time scale analysis: Determine the memory time. 

(vii) Second part of the time scale analysis: Take the Markovian limit and 
determine the typical time scale on which the relevant observables evolve; 
then, by comparing this time scale with the memory time, determine a 
posteriori whether the Markovian limit has been justified. In this manner, 
obtain physical criteria for the validity of the Markovian limit. 

(viii) Provided that the Markovian approximation is justified in the physical 
regime of interest, recover a Markovian equation of motion. If not, 

- solve a non-Markovian equation of motion; or 

- systematically expand around the Markovian limit; or 

- extend the level of description, then repeat the time scale analysis. 

Most quantum transport equations, notably rate equations and the 
quantum Boltzmann equation, refer to the evolution of single-particle 
distributions. Relevant are thus certain number operators 



Associated with this single-particle level of description is a relevant part of 
the statistical operator of the form 



the Lagrange parameters being adjusted so as to yield, at each time t, the 
correct occupation numbers 



{N t } := {aja,-} 






m(t) := (Ni)(t) 
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It is commonly assumed that the states {\i}} constitute a complete orthonor- 
mal basis; hence the observable 



N:=Y,Ni (87) 



represents the total number of particles. Furthermore, it is assumed that there 
are "no initial correlations," i. e., that 

P(t ) = Prel(t ) ■ (88) 



The resulting transport equation has the structure 

hi(t) = F[{ nj (t')\t < t' < t}] , (89) 

which only in the Markovian limit simplifies to 

h ( r\t) = F^[{n J (t)}} • (90) 



Such a Markovian limit is often tacitly assumed without rigorous justification. 
Here we wish to proceed more carefully: we first employ the projection method 
to derive the exact non-Markovian equation of motion; then we perform a 
thorough time scale analysis; obtain criteria for the validity of the Markovian 
limit; and only where appropriate, recover the usual Markovian expression. 

Different quantum transport equations, or individual terms therein, differ from 
each other in the choice of basis states {\i}} or in the form of the underlying 
microscopic dynamics. In the following sections we will discuss four specific 
microscopic processes -acceleration, level transitions, spontaneous pair cre- 
ation and binary collisions- and investigate how they manifest themselves in 
quantum transport equations. The examples are summarized in Table 2. 



For later reference we collect some useful formulae. The field opera- 
tors satisfy commutation or anticommutation relations, 

[a^a^ = [a],a]] T = , (91) 
[a,i,a]] T = 5ij , (92) 

depending on whether they describe bosons (upper sign) or fermions (lower 
sign). Since the relevant part of the statistical operator has the so-called 
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Table 2 

Four examples of quantum transport equations or parts thereof. The table lists 
the underlying microscopic dynamics, the choice of basis states (|p)=momentum 
eigenstates, |eQ)=eigenstates of the unperturbed Hamiltonian), the choice of the 
projector (M=Mori, KG=Kawasaki-Gunton), the interesting term in the general 
equation of motion (f=first term, m=memory term), and the resulting transport 
equation or part thereof (QBE=quantum Boltzmann equation). 



Microscopic dynamics 


10 


V 


term 


gives rise to 


acceleration in 
electric field 


|P> 




f 


acceleration term in QBE 


level transitions due to 
perturbing potential 


K) 


M/KG 


m 


rate equation for 
elastic scattering 


spontaneous pair creation 
in electric field 


Ip> 


M/KG 


m 


source term in QBE 


binary collisions 


Ip> 


KG 


m 


collision term in QBE 



Hartree form (cf. table 1), arbitrary expectation values of products of field 
operators can be calculated with the help of Wick's theorem [67]: 



(ABC...F) rel = [A'B'C" ...F'"} + [A* B" C\ . . F'"} + ... . (93) 



Here we have employed several shorthand notations. A, B, . . . represent arbi- 
trary annihilation or creation operators; and the "bullets" • indicate contrac- 
tions, defined by 

A'B* := (AB) re , . (94) 

Furthermore, it is understood that, e.g., 

[A•£••C•...F•••] = ±[A*C*B" ...F"*} , (95) 



the sign depending on whether the particles under consideration are bosons 
or fermions. More generally, each rearrangement r of the field operators must 
be accompanied by a factor 

e(r) := (±l) de ^ (96) 



which is determined by the degree of the permutation involved. Note that not 
all permutations are permitted: the order of two operators which belong to the 
same contraction must be maintained. Wick's theorem can then be written in 
the symbolic form 

(TT • • -)rei = g( r ) TT contractions . (97) 
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Each individual contraction may be evaluated using 

(ajaj)rei = Sijrii = S tJ [e y =p l]" 1 . (98) 

We also note the correlation matrix 

C %] := (SNi; SN 3 ) rel = 8 %3 n % {\ ± m) . (99) 



J f .2 Acceleration Term 

Acceleration gives rise to a non-vanishing first term in the general 
equation of motion. Our first example is almost trivial. We consider the 
acceleration of charged particles in a homogeneous, time-independent electric 
field E. Relevant observables are the occupation numbers of momentum eigen- 
states {|p}}. When applied to a relevant number operator N(p), the evolution 
superoperator simply causes a shift in momentum, 

U(t,0)N(p) = N(p + qEr) , (100) 



with q denoting the electric charge of the particles. Relevant observables are 
thus being transformed into relevant observables; relevant and irrelevant de- 
grees of freedom are not being mixed. Therefore, independently of the choice 
of projector, the memory term vanishes because QCN(p) = 0. There remains 
only the first term 



h(p,t) = i(p(t)\V(t)CN(p)) 

d (p(t)\V(t)U(r,0)N(p)) 

T = 



dr 

A 

"ck 



n(p + qEr,t) 

T = 

-qE-V p n(p,t) . 



This is the well-known acceleration term in the quantum Boltzmann equa- 
tion. It represents an important example in which the first term of the general 
equation of motion does not vanish even though all relevant observables com- 
mute (cf. our remark in section 3.4.2.). Since there is no integration over past 
history involved, the acceleration term is always Markovian and a time scale 
analysis not required. 
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4-3 Rate Equation for Elastic Scattering 



4-3.1 Non-Markovian Equation of Motion 

The underlying microscopic processes are transitions between en- 
ergy levels due to a perturbing external potential. Relevant observ- 
ables are the occupation numbers of the unperturbed levels; the appropriate 
basis states {\i}} are thus given by the eigenstates {|cq}} of the unperturbed 
Hamiltonian. Since the perturbing potential is assumed to be external, and 
not the result of particle interactions, the full Hamiltonian has the form 

H = H^ + V = J2^N t + J2V kl al ai . (102) 

i kl 

(Our notation does not account for possible level degeneracies. However, a 
generalization of the formalism which includes such degeneracies is straightfor- 
ward.) As the total number of particles is a constant of the motion, [H, N] = 0, 
any grand canonical state 

p eq [a, f3] := Z' 1 exp(-(3H + aN) (103) 
is stationary. 

In this particular application Kawasaki-Gunton and Mori projec- 
tor coincide. Observables which appear in this model are all of the single- 
particle form J20 tJ a]aj. It is therefore the effect on a\a 3 which is the only 
interesting property of a projector V . In this respect, the Kawasaki-Gunton 
projector "Pkg^'^)] anc ^ the Mori projector 7\[[a,0], defined in /? eq [a,0], are 
equivalent: both yield 

Va}a 3 =S t3 N t , (104) 

independent of or a. We may thus use a generic symbol V for the 

projector, and invoke its interpretation as either Kawasaki-Gunton or Mori 
projector only when it is needed or convenient. 

Evaluating the projected equation of motion yields a non-Markovian 
rate equation. Application of the interaction Liouvillian V to a field oper- 
ator a] yields (for both bosons and fermions) 

Vat = JT 1 5> w at , (105) 

k 
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and to a relevant observable Ni 



VN l = h- l Y J V kl (5uaia l -8 lk a\a l ) . (106) 



ki 



This implies VCNi = VVNi = 0, hence the first term in the equation of motion 
vanishes. So does the residual force term: Since the projector may be regarded 
as a Kawasaki-Gunton projector and since there are, by assumption, no initial 
correlations, the general argument applies (cf. section 3.4.2) according to which 
there is no residual force. Thus there remains only the memory term. Invoking 
now the interpretation of V as a Mori projector, we can take over the memory 
term of the Langevin-Mori equation: 

t-t 

».■(*) = £ / dT[-^(r)]Sn k (t-T) . (107) 



Here 

^{T) = [n(l±n)]- 1 (VN k ;T(0,T)VN l ) a , (108) 
Sn k (t) = n k (t)-n (109) 

and 

n:=(N k ) a VA; , (110) 



all expectation values and canonical correlation functions being evaluated in 
the stationary state /? eq [a,0]. Using 

£7?(r) = , (111) 

k 



the equation of motion may also be cast into the form 

t-t 

h i (t) = Y, J dr[- 7 f(r)]-K(t-r)-n 8 (t-r)] , (112) 

k n 



a form which already looks very much like a rate equation. It is, however, still 
non- M ar kovi an . 
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The memory matrix is calculated in perturbation theory. In case the 
perturbing potential is weak, the memory matrix may be evaluated to lowest 
non-trivial (i. e., second) order: 

, (2/h 2 ) Re [E« 8 \V u \ 2 exp(iu H T)} : k = i 
7f (r) = 2 L ¥ h • (H3) 

-(2/h 2 )Re [\V kt \ 2 exp(iu kt T)\ : k ^ i 



(Here we used the abbreviation ujki '■= ( e o — e o)/^0 H 1S lhi s resui t from second 
order perturbation theory which we want to subject to a time scale analysis. 



J f .3.2 Memory Time 



The memory time is related to the inverse coupling width. For our 

analysis we assume that the number of levels is large; that the system is close 
to its equilibrium configuration {n^ q }; and, more stringently, that only one 
particular occupation number rij(t) is perturbed appreciably from its equilib- 
rium value: 



+ Srij(t) : i = j 

eq ■ I ■ 

n l ■■ 1 T 3 



dt) = { ' n/ ' ■ (H4) 



The rate equation then describes the relaxation of rij(t) towards equilibrium: 

t-t 

Sn 3 (t) = - [ dT^(r)Sn 3 (t-T) . (115) 



From this expression one can extract both the memory time T^ em and the re- 
laxation time r^ el which are associated with the occupation number n r (Note 
that for each occupation number there is a different set of time scales.) Infor- 
mation pertaining to these scales is encoded entirely in the memory matrix 
element 7j(t). Inspection of the perturbative result shows that 7j(t) is essen- 
tially the Fourier transform of |V' J '| 2 with respect to the frequency difference 
u>ij. Hence if, by virtue of the perturbing potential, the level j couples to an 
entire "band" of other levels {/}, the band having a typical frequency width 
Aj, then 7j(i") is some distribution in r with characteristic width 1/Aj. The 
latter immediately determines the memory time: 

r^-l/A, . (116) 



As an example we consider electron transport in a conductor at low temper- 
atures (cf. section 4.3.5). The electrons move with Fermi velocity vp under 
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the influence of randomly distributed impurities. Each impurity has a typi- 
cal size r/ } which we assume to be large compared to the Fermi wavelength. 
The associated potential then couples states within a frequency band of width 
A ~ A(p 2 /2mh) ~ vpAk ~ vpjr\. The memory time is thus r mem ~ rj/vp, 
i. e., the typical time needed for an electron to pass through the interaction 
range of one impurity [35]. 



J f .3.3 Markovian Limit and Relaxation Time 

The relaxation time is related to the total transition rate. Provided 
the memory time is sufficiently small, one may take the Markovian and qua- 
sistationary limits. The general rate equation then becomes 

nt{t) = ^r kl [n k {t) - n^t)] , (117) 



where we have identified the transition rates 

00 

r kl := |dr[- 7 f(r)] , k±i . (118) 


Inserting the perturbative result for the memory matrix yields 

r kl = (27T/h 2 )8(u kl )\V ki \ 2 , k^i , (119) 

in agreement with Fermi's golden rule [68]. In the Markovian limit the relax- 
ation towards equilibrium is exponential, 

S' nj (t) = -rf'Sn^t) , (120) 
with the decay constant furnished by the total transition rate 



r tot 

3 



£ny • (121) 



From this we identify the relaxation time 

i = l/r!- ot • (122) 



The Markovian limit is justified for narrow resonances. With hind- 
sight we may now conclude that the Markovian limit is justified whenever 
the relaxation time is much larger than the memory time; or, equivalently, 
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whenever the total transition rate is much smaller than the coupling width 
("narrow resonances") [69,70]. While the memory time does not depend on 
the absolute strength of the potential, but only on its shape, the relaxation 
time increases as the potential becomes weaker. Therefore, the quality of the 
Markovian approximation -like the quality of perturbation theory- generally 
improves as \V kl \ — > 0. It is important to keep in mind, however, that despite 
this common feature the Markovian limit and perturbation theory represent 
two independent approximations. 



4.3.4 Dynamical Correlations 

Dynamical correlations do not depend on the state of the system. 

From our previous considerations we know that in the Markovian limit 

r tot . _ ^ 

7<,-(")H 1 Wu ■ ( 123 ) 

-r ki : k ^ 1 



Provided the total transition rate is much larger than individual rates, the 
memory matrix is almost diagonal; we may then set 

7J,-("Wk ■ lAii Vcu . (124) 



Inserting this approximation into Eq. (77) yields the Fourier transformed dy- 
namical correlation matrix: 

"k /, \ ^ r 'rel 



— ^-r . (125) 



ILOT 



rel 



This result depends only on relaxation time and frequency, but not on the 
thermodynamic state (temperature, density, etc.) of the system. 



4-3.5 Application: Drude Formula for the Electrical Conductivity 

The conductivity is defined in linear response theory. Macroscopic 
(normal) conductors at low temperatures owe their non-zero resistance to the 
presence of impurities [35]. According to the Kubo formula of linear response 
theory, their conductance G is determined by the current-current (/-/) corre- 
lation function: 

G(u) = (3(I;g < (u)I) eq . (126) 
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Here £J< is the Green's function superoperator (which includes the effects of 
elastic impurity scattering), and (3 denotes the inverse temperature. Instead 
of the global conductance one may also consider the local conductivity tensor 
CFik } which is essentially given by the autocorrelation function of the current 
density [60]. The current density, in turn, is related to the particle momenta 
such that 

a lk {u) = ^-{8P l] g < {u)8P k ) eOL , i,k = 1...3 . (127) 



Here i,k label the components of the conductivity tensor; e denotes the elec- 
tron charge, m the electron mass, V the volume; and SP is defined as 

£P:=P-(P) eq , (128) 



with P denoting the operator of total momentum. 



The conductivity is related to dynamical correlations of occupation 
numbers. At this point it seems very suggestive to identify the {Pi} as 
"relevant observables," and to evaluate their dynamical correlation function 
directly with the help of Langevin-Mori theory (Eq. (77)). However, the level 
of description spanned by these few observables is generally not Markovian; 
and as a result, the evaluation of the memory matrix may pose difficult prob- 
lems. In order to make use of the powerful Markovian approximation, one 
should rather extend the level of description. Provided the impurity density 
is sufficiently small and hence the perturbing potential sufficiently weak, we 
know from our previous considerations that elastic scattering of electrons off 
these impurities may be described by a Markovian rate equation. In this rate 
equation the place of the unperturbed levels | e ) is taken by the momentum 
eigenstates |p), and that of the occupation numbers {A;} by {A^(p)}. It is 
therefore the number operators {A^(p)} which span the appropriate level of 
description: this level is Markovian, and since 

^• = E^(P) > (129) 
p 



it contains the original level of description as a subspace. Using (75), the con- 
ductivity tensor can then be related to the dynamical correlation of occupation 
numbers: 

e 2 (3 

Vik(u) = — J{7Z)Z)Z)*PfcCqq'"<pH • ( 13 °) 
m V D a n' 
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Explicit evaluation yields the Drude formula. Since both the ordinary 
correlation matrix C and the dynamical correlation matrix S are diagonal, two 
of the three sums immediately collapse. Furthermore, at T — > 0, the diagonal 
elements of C are given by partial derivatives of the occupation numbers with 
respect to the Fermi energy: 



These derivatives are non-zero only on the Fermi surface. Hence -provided the 
impurity scattering is isotropic- we may replace 



where r no longer depends on the momentum p and denotes the elastic scat- 
tering time of electrons on the Fermi surface. The remaining summation over 
p is then performed easily: Using 



where n denotes the number density of electrons, we obtain the well-known 
Drude result [71-73] 






(132) 




(133) 




(134) 
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4-4 Source Term 



4-4-1 Non-Markovian Equation of Motion 

Our investigation is based on a simple model from quantum electro- 
dynamics. Strong external fields not only cause the acceleration of charged 
particles, but in addition may lead to the spontaneous creation of particle- 
antiparticle pairs. As an example we consider the spontaneous creation of 
e + e~-pairs in a homogeneous, time-independent electric field E, a process of- 
ten referred to as the Schwinger mechanism [74-78]. Relevant observables are 
the number operators 

N_(p,m z ) :=a\p,m z )a(p,m z ) , 

N + (p,m z ) :=b\p,m z )b(p,m z ) , (135) 
whose expectation values 

n T (p,m z ,t) := (N T (p,m z ))(t) (136) 



describe the momentum and spin distribution of electrons and positrons, re- 
spectively. Here a and b denote electron and positron field operators, p the 
momentum and m z the spin component. The transport equation will have the 
structure 

h T (p,m z ,t) ± qE ■ V p n T (p,m z ,t) = n s ° u (p,m,,t) , (137) 



with the known acceleration term (cf. section 4.2) and a new source term 
which accounts for spontaneous pair creation. In the following we will show 
how this source term can be derived from first principles. 

Our starting point is the Dirac equation 

ih±\tl>(t)) = H\tl>(t)) (138) 



with the Hamiltonian 

H = p ■ a + mfj + qA (139) 

and electrostatic potential 

A (r) = -E • r (140) 
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(where q = — \e\ for electrons). For later use we introduce some definitions: 
the time-dependent momentum 



p(t) := p + qEt 



:uv, 



the transverse energy 



(142) 



the total kinetic energy 



(143) 



and the dynamical phase 



1 

4> ft :=-j-jdt'e[p(f)] . 



(144) 



In our notation "longitudinal" 
direction of the electric field. 

The basis states 

\i,±) = \p(ti),m z ,±) , 



and "transverse" (_L) always refer to the 



(145) 



which correspond to momentum p spin component m z and positive or 
negative energy ±e[p(t 8 )], evolve according to [79-82] 



U(t f ,U 



t 



V 



t 



( e +t ^f* 



\-P*fi a *fi V e " % 



l/,+> 

I/,") 



(146) 



The evolution thus mixes positive and negative energy states, with respective 
amplitudes aji and These amplitudes are determined by the differential 
equation [83] 



<Xfi 



Pfi 



qE e_i_ 



- e - t2 ^f 

e +^f* 



/ 



hi 



(147) 



with initial conditions an = 1 and fin = 0, and the overdot indicating differ- 
entiation with respect totj. 
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The microscopic dynamics is described by a time-dependent Bo- 
goliubov transformation. In view of applying the projection method, the 
above results have to be translated into the language of field operators. Em- 
ploying the shorthand notation 



a(i,±) = a(p(t l ),m z ,±) 



:i48) 



for particle field operators, and making use of the general rule 

U(t h t i )a\il>) = a\U(t h t i )il>) , (149) 



one finds 




t A 



a f (/,+) 

« f (/," 



; (150) 



the evolution law for a follows by Hermitian conjugation. Identifying negative 
energy particle operators with positive energy antiparticle operators makes 
evident that pair creation is described by a time-dependent Bogoliubov trans- 
formation [67,84,85]. The modulus squared of the associated amplitude, |/3/;| 2 , 
equals the probability for having created an e + e~-pair with final momenta 
±p(t/) during the time interval 



From the evolution superoperator one obtains the Liouvillian 



r ■ 9 
dt f 



tt=t. 



which may be written as the sum 

c = c acc + V 



(152) 



of a diagonal part £ a cc, responsible for acceleration, and an off-diagonal part 
V which is responsible for the mixture of particle and antiparticle states and 
hence for pair creation. With the definition 



t f =t, 



(153) 



the latter is given by 



V 



a T U,+ 




a T U,+ 



(154) 
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The source term represents a particular example of a non-Markovian 
rate equation. Comparing Eq. (154) with Eq. (105), we note that with the 
identification 

ifci ^ 7T V'- )(, '' +) (155) 

there is a one-to-one correspondence between the microscopic dynamics of 
spontaneous pair creation and level transitions in an external potential. From 
this we infer that the source term must have the structure of a -generally 
non-Markovian- rate equation. Provided the electric field is sufficiently weak, 
E <C m 2 /hq } then \(3\ <C \(f>\ and hence V may be regarded as a small per- 
turbation. In this case we can immediately take over both Eq. (112) and the 
perturbative result (113); with the sole modification, however, that the evolu- 
tion superoperator 

T(0,r) exp(iQC acc Qr) (156) 

contained in the memory term now generates not only a phase factor [exp^o^r)], 
but also a shift in momentum. Therefore, the correct rate equation reads 

t-t 

n sou (p,m„±;t) = 2Re J dr /T(-r, -r)e- 82 ^- T '°)/?(0, 0) 

o 

X KP - <?Et, m z , =F; t - r) - n(p - qEr, m z , ±; t - r)] 

(157) 

with (3(tj,ti) = fiji and c/)(t/,t 8 ) = (f>ji. Upon identifying electron and positron 
occupation numbers 

n_(p,m z ,t) = n(p,m z ,+;t) , 

n + (p,m z ,t) = 1 - n(-p, -m z , -; t) (158) 
we obtain the source term [45] 

t-t 

n s q= ou (±p,±m„t) = 2Re J dr /3*(-r, -r) e - 82 ^- T '°)/3(0, 0) 

o 

x [1 — n_ (p — gEr, m z , t — r) 

-n + (-p + qEr, -m z ,t - r)] . (159) 

This source term involves an integration over the entire history of the system, 
thus accounting for finite memory effects and rendering the evolution of the 
occupation numbers generally non-Markovian. It exhibits two characteristic 
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time scales: (i) the memory time r mem (p), which corresponds to the temporal 
extent of each individual creation process and which indicates how far back 
into the past one has to reach in order to predict future occupation numbers; 
and (ii) the production interval r proc i(p), the inverse of the production rate, 
which corresponds to the average time that elapses between creation processes 
and thus constitutes the typical time scale on which the occupation numbers 
evolve. Both time scales refer to the specific momentum p with which the 
particles are produced. Only if r mem <C r proc i can memory effects be neglected 
and the evolution be considered approximately Markovian. 



The source term is consistent with the Schwinger formula. For sim- 
plicity, we will assume that (i) the field is weak 1), (ii) the system is 
dilute (n T 0), and (iii) t — > — oo. In the weak field limit the differential 
equation (147) implies 

fj ff exp(i2^ fi ) =fj ft (160) 

which, together with the decomposition 

<K-r,0) = <K-r,-oo)-<K0,-oo) > (161) 



reduces the source term to 



p,m,) = 2Re / dr /T(-r, -oo)/?(0, -oo) 



2Re |/T(0,-oo)/?(0,-oo 
_d_ 

d7 



T, -CO 



(162) 



T = 



The amplitude (3(tj,ti) depends on the momentum p through the initial con- 
dition p(t = 0) = p. In order to obtain the total electron production rate per 
unit volume, we have to consider 



w:=(2/h 6 ) J d>n s _ ou (p,m^) 

OO oo 

f d(el) ( d P ||n s _ ou (p,m,) . (163) 



The numerical factor 2 corresponds to the two possible spin orientations. 
Noticing that in the source term 
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_d_ 



r, -oo 



qE 



T = 



d 

dp||(l 
d 



P||( T )=P|| 



^[PllWjPlK- 00 )] 



(164) 



the integration over the longitudinal momentum becomes straightforward: 



dp\\ n s _ ou (p, m 2 ) = qE ■ |/?(oo, -oo)| 5 



(165) 



The modulus squared, |/3(oo, — oo)| 2 , represents the probability that for a given 
transverse momentum p_i_ a non-adiabatic transition between the energy levels 
±e[p(t)] will ever take place. It is given by the Landau-Zener formula [86-89] 



ire i 



oo, — oo 



exp 



kqE 



(166) 



The subsequent integration over (e^_) is then easily performed, yielding the 
final result 



(qE) 2 I irm 
w = — t-tt exp 



2\ 



4:lT 3 fl 



kqE J 



(167) 



This does indeed agree with the leading term in the Schwinger formula [76]. 



4-4-2 Memory Time 

The memory time combines classical and quantum mechanical time 
scales. In the weak field limit the memory time can be easily extracted from 
the source term (159). The factor 



-r, —t oc 



(eJqE) 



P \\/qEy + (eJqEf 



:i68) 



constitutes a Lorentz distribution in r, centered around p\\/qE with width 
e±/qE. Significant contributions to the source term thus come from times r 
which are smaller than (p\\ + e±)/qE. The typical magnitude of p\\ may be 
inferred if one views pair creation as a tunneling process from the negative 
to the positive energy continuum [90] (cf. Fig. 2): the barrier between these 
continua has a spatial width of the order e±/qE } inducing a corresponding 
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Fig. 2. The spectrum of the free Dirac-Hamiltonian is comprised of a positive and a 
negative energy continuum, separated by a gap of magnitude 2m. (This gap widens 
to 2e_|_ if one fixes the transverse momentum.) In the ground state the negative 
energy continuum ("Dirac sea") is completely filled. Applying a homogeneous elec- 
tric field "skews" these energy levels, due to the linearly rising potential. It is then 
possible for an electron to tunnel from the lower to the upper continuum, a process 
which corresponds to pair creation. 

momentum scale Ap|| ~ hqE/e±. With this typical scale one obtains 



The memory time combines two time scales of different origin: (i) The time 
h/c±_ is proportional to h and therefore of quantum mechanical origin. It cor- 
responds -via the time-energy uncertainty relation- to the time needed to 
create a virtual particle-antiparticle pair, and may thus be regarded as the 
"time between two production attempts." (ii) The time e±/qE } on the other 
hand, is independent of h and therefore classical. It can be interpreted in var- 
ious ways, depending on the picture employed to visualize the pair creation 
process. If pair creation is viewed as a tunneling process (Fig. 2), the classical 
memory time coincides with the time needed for the wave function to traverse 
the barrier with the speed of light [91]. Alternatively, pair creation may be 
viewed as a non-adiabatic transition between the two time-dependent energy 
levels ±e[p(t)] (Fig. 3). In that case the classical memory time corresponds to 
the width of the transition region, i.e., the region of closest approach of the 
two levels. 

For weak fields, E <C m 2 /hq } the classical memory time dominates: 



mem 



k £_|_ 

e_i_ qE 



(169) 




e_L h 



(170) 
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lenergy 



asymptote: 




asymptote: 
\s-qE x 



Fig. 3. Pair creation viewed as a non-adiabatic transition between two 
time-dependent energy levels. The temporal extent of the transition region is of 
the order ej_/qE. 

4.J f .3 Markovian Limit and Production Interval 

The production interval displays a characteristic E~ l exp(c/£')-de- 
pendence. In the weak field limit we can give a simple estimate for the 
order of magnitude of the production rate. According to Eq. (162) 



n sou ~0(|/3/3|) 



:i7h 



which together with the differential equation (147), 



~ 0(qE/e ± ) 



(172) 



and the Landau-Zener formula (166), 



exp 



7T£ 1 



2nqE J 



(173) 



yields the estimate 



• sou 

n ~ 



ire 1 



2nqE 



(174) 
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Inversion gives the typical production interval 




) 



(175) 



Pair creation in strong fields is always non-Markovian. The evolu- 
tion of occupation numbers is Markovian whenever the production interval far 
exceeds the memory time, r proc i ^> r mem . This is the case only as long as the 
external field is much weaker than the critical value m 2 /hq. As soon as the 
field strength reaches or exceeds that critical value, the situation changes: The 
source term (159), which was derived in the weak-field limit, is then merely 
a rough estimate. Already this weak-field estimate becomes non-Markovian, 
indicating that at this point any conventional -i. e., Markovian- transport 
theory must break down. (The criteria for the validity of second order per- 
turbation theory and of the Markovian limit are thus not independent: both 
approximations break down at the same critical field strength. This coinci- 
dence might have been expected, given that the model features only a single 
adjustable parameter (E).) In other words, our analysis has revealed that 
in the strong-field domain pair creation must always be described by a non- 
Markovian transport theory. This non-Markovian nature may have important 
ramifications. For example, it may lead to oscillations of the relevant entropy, 
and thus temporary violations of the _ff-theorem, on the scale of the memory 
time (cf. appendix and Refs. [45,92]). 

4-5 Collision Term 

4-5.1 Non-Markovian Equation of Motion 

The underlying microscopic processes are binary collisions. We con- 
sider the evolution of a system of identical particles (bosons or fermions) which 
interact through binary (2 — > 2) collisions. Relevant observables are the oc- 
cupation numbers of momentum eigenstates {|p)}. For brevity, we will often 
write 



We assume the system to be confined to a large, but finite, volume so that 
we can use the box normalization 




(176) 




(177) 
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and completeness property 

£ I0<i| = l • (178) 



Our notation does not account for spin or other discrete quantum numbers. 
However, a generalization of the formalism which includes such additional 
quantum numbers is straightforward. 

The microscopic Hamiltonian has the form 

H = + V = J2 <LNi + i p/fc| V| jO±«Mw • (179) 

i ijkl 



Here e^ in denotes the kinetic energy associated with momentum p ? , and 

(lk\V\ji)± :=(lk\V\ji)±(kl\V\ji) (180) 



denotes an (anti-)symmetrized matrix element, the sign depending on whether 
the particles are bosons (upper sign) or fermions (lower sign). Defining 

Vik\ji := {lk\V\ji)±a\a\a 3 a % (181) 



with the symmetries 

= V k l\ji = \l k \ij , \ lk y Jt = Vij\ k l , 



VjMii — Vkttii — VjMii , V,u ■■ — Vi;.\kl , (182) 



one can also write 



V = i^V m . (183) 

ijkl 



We assume that the dynamics is invariant under spatial displacements, which 
implies momentum conservation: 

Vik\ji = Spi+pjVk+ptViktfi ■ (184) 



For future reference we note that 

[Vik\ji, N{p)} = {Sp tP + 8 PjP - 8 PkP - 8 PlP )V m . (185) 
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This commutator vanishes whenever Vik\ji corresponds to zero momentum 
transfer. For example, if p; = p 3 and hence p/; = p 8 due to momentum 
conservation, then the Kronecker symbols add up to zero. More generally, 

[V lkl3t , N(p)] = if I = j, I = i, k = j, or k = i . (186) 



We employ the Kawasaki-Gunton projector. In this application the 
Kawasaki-Gunton and Mori projectors do not coincide. The choice of the 
Kawasaki-Gunton projector is motivated by its three potential advantages 
discussed in section 3.4.2: vanishing residual force, vanishing first term, and 
-perhaps most importantly- applicability of Wick's theorem. Thus our start- 
ing point will be the Robertson equation (63). That the first term vanishes, 
should be checked explicitly: Application of the interaction Liouvillian to a 
relevant observable N(p) yields 

VN(p) = (4h)- 1 J2lVm,N(p)] . (187) 

ijkl 

According to Wick's theorem, only those terms in the sum can have a non- 
zero expectation value (. . .) re i for which {/,&} = {j, i}; yet in exactly that 
case, according to Eq. (186), the commutator vanishes. Hence no contribution 
survives, the first term must indeed be zero: 

i{£N(p)) rem = i{VN(p)) rem = . (188) 

This in turn implies 

V{t)CV{t) = V(t)VV(t) = . (189) 

In the Robertson equation, therefore, there remains only the memory term; 
and in this memory term all complementary projectors (Q) may be omitted. 



The Robertson equation gives rise to a non-Markovian collision 
term. We will evaluate the Robertson equation in second order perturbation 
theory. Combining the Robertson equation (63) with the general perturbative 
result (51) yields, to second order, 

t-t 

n{p,t) = - j dr(V^°)(0,r)ViV(p)) rel(t _ T ) . (190) 
o 
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Using now Eqs. (187) and (185), as well as the property 

W (o) (0,t)V^| J8 = expiiSujik^Vik^ (191) 



with 

r Hk\-ji 
leads to 



8u m := % '(eLin + 4in - 4in ~ 4in) > ( 192 ) 



q abed ijkl 

><([Vd c |6 a , VJjt|ji]) re l(t_ T ) • (193) 



The expectation value in the integrand can be calculated with the help of 
Wick's theorem (a tedious, but straightforward calculation). One finds 

E([^c|6a, V m ]} = 4|</£;|V|jz) ± | 2 [n 8 n J (l ± n k )(l ± n,) - (ij ^ A;/)] (194) 

abed 



where n 8 - := (aja 4 -), all expectation values being evaluated in the state p Te \(t — r) 
Of the remaining four summations (over i } j, A;, /) one more can be performed, 
yielding finally the non-Markovian collision term 

t-t 

h(p,t) = h- 2 Re f dr /[PP3|P2Pi](r)-C[pp 3 |p2Pi](<-r) (195) 

n P1P2P3 



with collision bracket 

C[PP3|p 2 Pi](0 : = n(p 1 ,t')n(p 2 ,t / )[l±n(p3,t / )][l±n(p,t / )] 

-(pip 2 ^p 3 p) (196) 

and prefactor 

/[PP3|p2Pi](r) :=exp(z^[pp 3 |p 2 pi]r) • |(pp 3 | UIP2P1 )± | 2 £ p+P 3,P2+Pi . 

(197) 

Here we have inserted an extra Kronecker symbol to make the conservation of 
momentum explicit. The frequency increment 8u which appears in the phase 
factor is essentially the net "energy gain" in a collision 1 + 2 — > 3 + X. Of 
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course, it must vanish if the collision conserves energy. Energy conservation, 
however, has not yet been established. Therefore the phase factor exp(i8u . . .) 
may not be neglected; in fact, it will play a crucial role in the analysis of time 
scales. 



J f .5.2 Memory Time 

For our analysis we make some simplifying assumptions, (i) The 

dynamics is nonrelativistic; hence for particles with mass m, 

M PP3 |P2Pi] = ^t(p 2 + P3-P2-pD • (198) 

Zmn 



Defining the momentum transfer 

hq : = p - pi = p 2 - p 3 (199) 

and the velocity change 

v := (p - p 2 )/m = (pi - p 3 )/ra , (200) 

one may also write 

Su = q-v . (201) 



(ii) We assume that the interaction has a typical range ra and make the generic 
choice of a Gaussian potential [38] 

V(r) oc exp(-r 2 /ri 2 ) . (202) 



Neglecting the effects of (anti-)symmetrization of the matrix elements, this 
leads to 

KPP3|V|p 2 p 1 ) ± | 2 ocexp(- (? V/2) • (203) 



(iii) The system is dilute so that one may neglect the enhancement /blocking 
factors [lira], (iv) The system is close to its equilibrium configuration {ra eq (p)}, 
i.e., 

ra(p,t) = ra eq (p) + Sn(p,t) (204) 
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where the deviations 8n are small. More stringently, we assume that only one 
particular occupation number n(p , t) is perturbed from its equilibrium value: 

f ^ if p = p 

8n(p,t)l T . (205) 

I = otherwise 

In the collision term only the deviation of the collision bracket from equilib- 
rium, SC := C — C eq , contributes. For p = p , 

<fC[poP 3 |p2Pi](0 = -n eq (p 3 )^(po,t / ) . (206) 

(Variations of ra(pi) and n(p 2 ) do not contribute because if, for example, pi = 
p then p 2 = P3 due to momentum conservation and C = C[poP3|p3Po] = 0.) 
The collision term then describes the relaxation of n(p ,t) towards equilib- 
rium: 



t-t 

8n(p 0} t) = -h~ 2 Re / dr ^ exp(i8uT)\{p p 3 \V\p 2 pi}±\ 2 

P1P2P3 

x^Po+P3,Pi+p 2 n e(i (p 3 )8n(p ,t - r) . (207) 

(v) Finally we make an assumption about the shape of the equilibrium dis- 
tribution. Again, we neglect the effect of (anti-)symmetrization and choose a 
Maxwell-Boltzmann distribution with inverse temperature (3: 

n eq (p 3 ) oc exp(-(3p 2 3 /2m) . (208) 



The memory time corresponds to the typical duration of one indi- 
vidual collision. From the relaxation equation (207) one can extract both 
the memory time r mem (p ) and the relaxation time r re i(p ), associated with 
the occupation number n(p ,t). Strictly speaking, there is a different set of 
time scales {r mem (p), r re i(p)} for each observable n(p,t). We will be content, 
however, with performing the analysis in the special case p = 0. Other mem- 
ory and relaxation times can be found with the same procedure. In fact, it 
can be shown -but will not be shown here- that any r mem (p) can never be 
greater than r mem (0) [38]. Hence r mem (0) is actually an upper limit for the 
various memory times, and it is fair to regard r mem (0) as the memory time of 
the system: r mem (0) = r mem . Likewise, the relaxation time r re i(0) may be re- 
garded as representative. And although by definition r mem and r re i only apply 
to small deviations from equilibrium, they are nevertheless useful measures for 
the typical time scales on which the relevant and irrelevant degrees of freedom 
evolve, even in situations further removed from equilibrium. 
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In order to evaluate the relaxation equation explicitly, we (i) assume that the 
volume is large, permitting to take the continuum limit 

E ^ ^ / ^ , <W> 7^ 3 (p.- -Pi) ; ( 209 ) 

(ii) absorb the ^-function into the integration over p 3 ; (iii) substitute (pi, p 2 ) — > 
(q, v); and (iv) use p = p = 0. Up to a normalization factor we obtain 



t-t 



8n(0, t) oc Re / dr / d 3 v / d 3 q exp 



1 2 2 

zq • vr — — q i] - 



„ ( h 2 q 2 mV 
\ 2m 2 i 



8n(0,t-T) . 



(210) 



The integration over v amounts to Fourier transforming a Gaussian of width 
(/3m) -1 / 2 , which yields another Gaussian distribution in (r + ih(3)q of width 
{fjm) l l 2 . One is left with 



t-t 



cm(0,t)ocRe /dr/d 3 gexp 



o 

t-t 



2 \ (3m m 



Sn(0,t- t) 



ocRe j <1t(t 2 + 2ih(3T + r] 2 (3m)- 3/2 8n(0,t - t) . (211) 
o 

Any normalization factors have been ignored, since we are only interested in 
the form of the integrand's r-dependence. Clearly, the integrand exhibits two 
characteristic time scales, (r]y/f3m) and (h(3). Times r which are greater than 
both of these characteristic scales yield only a negligible contribution to the 
integral. In other words, only those times r contribute significantly which are 
smaller than either (r]y/f3m) or (h(3). This finding may be summarized in a 
compact fashion by setting 



rj\l (3m + h(3 



(212) 



The two contributions to the memory time have very different origins: 

- The time r^ 1 ^ := rwj (3m is independent of % and therefore of entirely 
classical origin. Taking into account that (3 oc (m(u 2 )) -1 , one finds that 



class 

OC 



vN(v 2 ) (213) 
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is just the average time needed for a particle to pass through an interaction 
range. During that time correlations between colliding particles may not 
be neglected. The classical memory time depends both on the form of the 
interaction (through r/) and on the particle distribution (through (3). 
- The time T^" m := h(3, on the other hand, is proportional to h and therefore 
of purely quantum mechanical origin. It is obviously a remnant of the coher- 
ent superposition of individual 2 — > 2 scattering processes. More precisely, 
it describes the average temporal extent of a quantum mechanical scatter- 
ing process associated -via the time-energy uncertainty relation- with the 
characteristic energy scale fi~ l of the system. For shorter times r < h(3 en- 
ergy conservation is not yet established, the colliding particles are still off 
mass shell, and hence any subsequent collision would not be independent 
of the first; in short, for r < h(3 a scattering process cannot be regarded 
as completed. This quantum mechanical memory time ("off-shell time") is 
independent of the interaction and depends only on the particle distribution 
(through (3). 

In summary, the memory time corresponds to the typical duration -both clas- 
sical and quantum mechanical- of a single binary collision. As in the case of 
level transitions, the memory time does not depend on the strength, but only 
on the form (range r/) of the interaction. There is an additional feature, how- 
ever, in that the memory time now also depends (through (3) on the particle 
distribution. Both the classical and the quantum mechanical memory time 
decrease with increasing temperature. However, their respective temperature- 
dependence differs: \fj3 vs. (3. There is a "critical" temperature T c at which 
both memory times are of equal importance: 

kT c ~ h 2 /mri 2 . (214) 

For lower temperatures, T < T c , the quantum mechanical memory time dom- 
inates; while for higher temperatures, T > T c , the classical memory time 
prevails. Example: For a dilute gas of interacting nucleons (mc 2 ~ 1 GeV, 
j] ~ 1 fm) this "critical" temperature is kT c ~ 40 MeV. 

We finally mention that the above result for the memory time agrees with 
an analysis by Danielewicz who, instead of the projection method, employed 
Green's function techniques [38]. Earlier discussions of the collision term can 
be found in, e. g., Refs. [27,93]. 

J f .5.3 Markovian Limit and Relaxation Time 

The relaxation time corresponds to the average time that elapses 
between two successive collisions. Provided the memory time is suffi- 
ciently small, one can take the Markovian limit 8n(p 0} t — r) — > $n(p ,t), 
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as well as the quasistationary limit (t — to) — > oo. The r-integration in the 
relaxation equation (207) is then easily performed, yielding 

oo 

Re J dr exp(z<W) = ttS(Su) . (215) 
o 

The above approximations thus imply energy conservation in each individ- 
ual collision. According to Fermi's golden rule (119) the transition rate for 
scattering (p3,p) — > (Pi,P2) is given by 

r[pP3|p 2 Pi] = (2ir/h 2 )S(Su) |(pp 3 |V|p 2 pi }±| 2 £p +P 3,p 1+ p 2 • ( 216 ) 



Taking the sum over all final states (pi,p 2 ), dividing by the incident flux 
|v 3 — v|/0, and multiplying by an extra factor 1/2 to account for the identity 
of the colliding particles, yields the total cross section 

^tot[p,p 3 ] = i — - — i • \ J2 r [PPs|P2Pi] • (217) 
l v 3-v| 2 pip2 



With this expression for the total cross section the relaxation equation takes 
on the simple form 

Sn(p ,t) = ( l v 3 - v | cr tot [po,p 3 ] ^eq(ps) j Sn(p ,t) . (218) 



It has the structure of an exponential decay law. For p = v = 0, the decay 
rate is given by 



r = 75El v Vtot[0,p'] n eq (p') 

Si p' 



(219) 



where n := N/ft denotes the particle density and (...) the equilibrium average. 
Associated with this decay rate is the relaxation time 

r rel = — ^— . (220) 

n{ya) 



The relaxation time is intimately connected with the mean free path Xj = 
1/ncr; in fact, it is nothing but the average time between two successive colli- 
sions. 
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The Markovian limit is justified whenever the average time between 
collisions is much larger than the temporal extent of each individual 
scattering process. The Markovian limit becomes exact as V — >■ (and 
hence a — > 0) or n — > 0; i.e., for weakly interacting or very dilute systems. 
In contrast, the quality of the Markovian limit does not generally improve 
with increasing temperature, since the average velocity cancels from the ratio 
r mem/ r rei- It might improve nevertheless if -as is often the case- the cross 
section decreases with increasing temperature. 
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5 Conclusion 



In the preceding sections we have demonstrated the use of the projection 
method in several practical applications. Following the agenda laid out in 
section 4.1, we have employed the projection technique to investigate how a 
variety of physical processes, such as level transitions, binary collisions or spon- 
taneous pair creation, manifest themselves in macroscopic transport equations. 
In all these examples we have ascertained that the projection method permits 
the derivation of macroscopic transport equations from the underlying micro- 
scopic dynamics; that the thus obtained transport equations are generally not 
Markovian, but exhibit a finite memory; that these non-Markovian equations 
may be subjected to a rigorous time scale analysis, yielding criteria for the va- 
lidity of the Markovian limit; and that, where appropriate, the Markovian limit 
then allows one to recover conventional -i. e., Markovian- quantum transport 
equations. These findings substantiate the claims made in the introduction 
(section 1) about the advantages of the projection method. Furthermore, they 
confirm the general statement made in section 3.1: namely, that the mapping 
of the influence of irrelevant degrees of freedom onto a temporal non-locality, 
and the resultant possibility of analyzing and exploiting well-separated time 
scales, are the primary reasons for the success of the projection method. 

In concluding, we wish to point out some current problems in quantum trans- 
port theory which we consider worth investigating and which, as we believe, 
may lend themselves to a treatment with projection techniques: 

- the formulation of a transport theory for highly excited many-body systems 
that goes beyond the conventional Boltzmann description by accounting for 
memory effects, the existence of several particle species, and the possibility 
of cluster formation (i. e., non-trivial correlations); 

- the study of dissipative phenomena and of possible memory effects in effec- 
tive quantum field theories; 

- the study of the transport properties of quantum systems whose classical 
counterparts are chaotic; and 

- a better understanding of the phenomenon of decoherence. 

This list of potential projects is certainly influenced by our own research inter- 
ests and by no means complete. We hope that our review will help stimulate 
investigations both of the above and of many other related issues. 
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A Entropy and Irreversibility 



A.l Entropy and Relative Entropy 

The entropy characterizes the lack of information as to a system's 
microstate. Let A be a (classical, not quantum) random experiment with 
possible outcomes {a 8 } and respective probabilities {p(a 8 )}. An "ignorance 
measure" I[A] } which characterizes the lack of information as to the outcome 
of A, should have the following three properties: 

(i) I[A] is a continuous, symmetric function of the probabilities {p(a 8 )}. 

(ii) The lack of information is non-negative, I[A] > 0. Furthermore, I[A] 
vanishes only if the outcome of the experiment is certain, i. e., if {p(a 8 )} = 
{0,... ,0,1,0,...}. 

(iii) The lack of information is additive. Let B be a second experiment with 
possible outcomes {b 3 } } and A D B the combined experiment with out- 
comes {a 8 - fl b 3 }. Then the total lack of information, I[A D B], is the sum 
of I[A] and the expected lack of information as to B } given the outcome 
of A: 

I[AnB] = I[A] + '£p(a i )I[B\a i ] . (A.l) 

i 

Here the measure /[Afl-B] is defined in terms of jj9(a 8 n&j)}, I[A] in terms 
of {p(a 8 )}, and 7[_B|a 8 ] in terms of the conditional probabilities {p(6j|a 8 )}, 
the various probabilities being related via Bayes' rule 

p(a t n bj) = p(a t )p(b 1 \a t ) . (A. 2) 

According to a celebrated theorem by Shannon [94-96] the above properties 
specify the ignorance measure uniquely (up to a prefactor): 

I[A] = -k^2p(a l ) lnp(a l ) , k > . (A. 3) 

i 

Let {i} denote a collection of discrete classical microstates of a physical system 
and {pi} their respective probabilities. The physical entropy 

S[p}-= -k^pilnpi , (A.4) 

i 

with k chosen to be the Boltzmann constant, is then nothing but the amount 
of missing information as to the outcome of a most accurate experiment (i. e., 
an experiment which can resolve individual microstates). 



59 



The above definition of the entropy can be extended to quantum systems. Any 
statistical operator p may be written as an incoherent mixture of pure states, 

p = Y,Mi){i\ ■ (A.5) 



The lack of information as to the system's micro- (or pure) state is then given 
by S[p] } which in turn coincides with the von Neumann entropy 

S[p] ■= -ktr(plnp) . (A. 6) 



For continuous probability distributions the entropy is not well de- 
fined. For the entropy of a continuous (classical) probability distribution 
p(x) one could envision the following definition which, however, turns out to 
be inappropriate: One might discretize the distribution by dividing the sup- 
port manifold into bins {A;}, set 

Pi := J dx p(x) , (A. 7) 



and define 



vol(A 8 )^0 \ / 



Yet the limit is ill-defined: this "naive" entropy diverges. And even after sub- 
tracting the divergent contribution (— k ln[vol(A 4 -)]), the finite expression 

s Live\p] '■= ~ k / dxp(x)lnp(x) (A. 9) 



is still not acceptable, because it is not invariant under coordinate transfor- 
mations. A physical information measure should not depend on the particular 
coordinate system in which the distribution is represented. 



A more suitable concept is the relative entropy. Instead of asking 
"How much information as to the microstate is lacking?" one might also ask 
the reverse question: "How much information do we already have, relative to a 
state of total ignorance?" In the case of a discrete probability distribution the 
state of total ignorance corresponds to a uniform distribution {m 8 - = const.}. 
The information gained by replacing {m 8 } with a nonuniform distribution 
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{pi} is then given simply by the difference of entropies, (S'fm] — iSfp]). This 
difference may also be written in the form 



Sdi B cMp):=fc£p.- ln — ^° • ( A - 10 ) 

m; 



It is in the latter form that the information gain can be defined also for 
continuous distributions (with compact support): 

ScontMp) := k f dx p(x) In -^T > . (A.ll) 
J mix) 



Except for pathological cases, this integral is well-defined and finite. Further- 
more, it is invariant under coordinate transformations x — > y, provided not 
only p } but also m is transformed: 

/ dx p{x) In ^\ = [dy p(y) In (A.12) 
J mix) J m(v) 



where 

p(x)dx = p(y)dy , m(x)dx = m(y)dy . (A. 13) 

In arbitrary coordinate systems the state m of total ignorance no longer neces- 
sarily corresponds to a uniform distribution. The above definition of S con t(m\p) 
ensures that even in these coordinate systems the information gain is well- 
defined. 

The preceding considerations are intended to motivate the introduction of the 
concept of relative (or Kullback) entropy [97,98]. The relative entropy S(q\p) of 
a distribution p with respect to another distribution q is just defined as above 
-in both its discrete and continuous versions- but with q and p now represent- 
ing arbitrary distributions. This generalizes the concept of information gain 
to cases in which the distribution q no longer represents the state of total 
ignorance. The relative entropy is still non-negative, finite and coordinate- 
independent; it vanishes only if the two distributions coincide: 

S(q\p) = iff p = q . (A. 14) 



The definition of relative entropy can be readily extended to quantum systems. 
For two arbitrary statistical operators <r, p one defines 



S(a\p) := fctr^ln/o-liKT)] > . (A. 15) 
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Note that in general the relative entropy is neither symmetric, S(cr\p) ^ 
S(p\cr) } nor does it equal the difference of ordinary entropies, S(cr\p) ^ S[cr] — 
S[ P ]. 



The relative entropy measures the distance between distributions. 

To test the quality of a statistical hypothesis, one usually employs the so- 
called x 2 test. The quantity \ 2 characterizes how significantly the observed 
relative frequencies {pi} (data) deviate from the theoretical predictions 
(hypothesis); \ 2 may thus be regarded as a measure for the distance between 
these distributions. It is defined as 

X 2 :=N± { -^^ , (A.16) 



with r denoting the number of bins (or classes) and N the total number of 
data points (events). The distance measure \ 2 1S intimately connected with 
the relative entropy S(q\p). Provided the two distributions p } q are sufficiently 
close to each other, one may expand S(q\p) in powers of (pi — qi) and obtains, 
to lowest order [100], 

%W^4x 2 • (A-17) 



This result strongly suggests that the relative entropy, too, may be regarded 
as a measure for the distance between distributions. Note that except for 
infinitesimal distances the relative entropy, like % 2 , is generally not symmetric: 
S(q\p) ^ S(p\q). 



A. 2 Three Uses of Relative Entropy in Statistical Physics 



Entropy and relative entropy serve as diagnostic tools. Entropy and 
relative entropy are not physical observables like energy or particle number. 
In quantum mechanics there is nothing like a Hermitian "entropy operator." 
Rather, entropy and relative entropy serve primarily as mathematical tools to 
analyze and compare properties of different probability distributions. Below 
we will discuss three important uses of these tools: (i) finding unbiased distri- 
butions, (ii) testing the quality of approximations, and (iii) testing whether a 
level of description may be contracted. 
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A. 2.1 Finding Unbiased Prior Distributions 



Minimizing the information gain yields a generalized canonical dis- 
tribution. Often a probability distribution is not known exactly. Rather, 
only certain constraints are given of the form 



T,iPiG a (i) 

J dx p(x)G a (x) 

tT(pG a ) 



(G a ) 



(A.18) 



the particular form depending on whether the distribution is classical dis- 
crete, continuous, or quantum. From this limited information one likes to 
infer that distribution p which, while satisfying all constraints, is "maximally 
non-committal" with regard to missing information and hence "least biased." 
Mathematically, this requirement amounts to minimizing, under the given con- 
straints, the information gain S(m\p) relative to a state m of total ignorance 
[101-104]. 

In the discrete and quantum cases, minimizing S(m\p) is equivalent to maxi- 
mizing the ordinary entropy S[p\. In the continuous case, however, there is the 
additional complication of finding, in the given coordinates x, the -possibly 
nonuniform- distribution m(x) which represents total ignorance. This igno- 
rance distribution m(x) can usually be determined on the basis of symmetry 
considerations: If the structure of the support manifold is invariant under some 
characteristic group (e. g., translations, rotations, canonical transformations) 
then m(x)dx must coincide with the group-invariant measure [105,106]. In all 
cases minimizing the information gain yields a generalized canonical distribu- 
tion: 

p i = Z- 1 exp(-\ a G a (i)) , Z = £exp(-A a G a (0) 

i 

p(x) = Z~ l m(x) exp( — X a G a (x)) , Z = J dxm(x) exp( — X a G a (x)) 
p = Z- 1 exp(-A a G' a ) , Z = tr exp(-A a G a ) . (A. 19) 



(A summation over a is implied.) One obvious application of these generalized 
canonical distributions is to equilibrium statistical mechanics, where the {G a } 
represent the constants of the motion. 
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A. 2. 2 Testing the Quality of Approximations 



The relative entropy measures the size of the error. From now on 
we focus on quantum systems. Let the macrostate of a quantum system be 
characterized by the expectation values of the relevant observables {G a }; let 
{g a (t)} be the exact expectation values, and {g' a (t)} some approximate values 
obtained after, e. g., taking the Markovian limit or doing perturbation the- 
ory. Associated with both sets of expectation values are generalized canonical 
states p ca ,n[ga(t)] and p c &n[g' a {~t)\i respectively. An appropriate measure for the 
quality of the approximation g g' is furnished by the distance between 
these canonical states, i. e., by the relative entropy S(p ca _ n [g' a (t)]\p ca _ n [g a (t)]). 
Provided the deviations of the approximate from the exact expectation values, 

Ag a (t) := g' a (t) - g a (t) , (A.20) 

and of the approximate from the exact associated Lagrange parameters, 

A\ a (t) := \ a '(t) - \ a (t) , (A. 21) 

are sufficiently small, then the relative entropy is given by 

S(pUg' a (t)]\pU9a(t)]) = \C ab A\ a A\ h + 0((AXf) 

= \(C- l r h Ag a Ag b + 0({Agf) . (A.22) 

The relative entropy thus increases quadratically with the deviations AA or 
Ag. The "metric tensor" C is just the correlation matrix defined in Eq. (31). 
It may be evaluated (to lowest order in AA or Ag) either in the exact or in 
the approximate canonical state. 

A.2.3 Testing the Feasibility of Level Contractions 

Whether or not a level of description may be contracted, can be 
decided by comparing ordinary entropies. Sometimes the characteri- 
zation of the macrostate by the expectation values of {G a } is redundant: a 
smaller level of description, 

span{l, F a } C span{l, G a } , (A. 23) 

might already be sufficient, yielding approximately the same macrostate: 

Pc,n[fa(t)] ~ Pc,n[ga(t)] • (A. 24) 
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If this is the case then the level of description may be contracted from the 
original (span{l, G a }) to its subspace (span{l, F a }). An appropriate measure 
for the feasibility of such a contraction is again furnished by the relative en- 
tropy S(p ca , n [f a (t)]\p ca _ n [g a (t)]). One can show that when p ca , n [fa(t)] is obtained 
from p can [g a (t)] via level contraction, the relative entropy coincides with the 
difference of ordinary entropies: 



S{p ca n[fa(t)]\pc a n[g a (t)]) = S[p can [f a (t)]] - S[p can [g a (t)]] > . (A.25) 



Whether or not the level contraction span{l, G a } — > span{l, F a } is permitted, 
can thus be decided simply by comparing the respective ordinary entropies. 
Their difference should not exceed some prescribed upper bound which de- 
pends on the accuracy desired for the relevant expectation values (cf. the 



A. 3 Temporal Variation of Entropies 

Non-trivial statements refer to the relevant entropy. The von Neu- 
mann entropy associated with the full statistical operator of a closed quantum 
system, 



is constant in time. Hence for the study of irreversible phenomena it is useless. 
Non-trivial statements all refer to the time-dependent relevant entropy, 



which is the entropy associated with the relevant part of the statistical oper- 
ator. It measures, at all times t, the amount of missing information as to the 
pure state of the system if one is given the current expectation values {g a (t)} of 
only the relevant observables. The transition pit) — > p re \(t) } and consequently 
S — y S re \(t) } is often referred to as "coarse-graining." It always involves dis- 
carding information about irrelevant degrees of freedom, a truncation which 
is reflected in the general inequality 




S[p(t)]:=-ktT(p(t)lnp(t)) 



(A.26) 



Srel[g a (t)] ■= -ktr(p re l(t)llip re l(t)) 



(A.27) 



Sr e l[ga(t)]>S[p(t)] Vt 
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About the relevant entropy one can make two distinct, logically independent 
statements: 

(i) Second law. Whenever both the initial and the final macrostate of an iso- 
lated system are characterized by a generalized canonical distribution -in 
particular, whenever the evolution of the system leads from one equilib- 
rium state to another equilibrium state- then the associated final relevant 
entropy can never be smaller than the initial relevant entropy. This is true 
even if the initial and final levels of description differ. No statement is 
made about the behavior of the relevant entropy at intermediate times. 

(ii) H -theorem. If a level of description is Markovian then its accompany- 
ing relevant entropy increases monotonically at all times. (Boltzmann's 
original _ff-theorem is one particular manifestation of this more general 
result.) 

These two statements will be derived and discussed separately. 
A. 3.1 Second Law 

The final relevant entropy cannot be smaller than the initial relevant 
entropy. Let {G a } and be two -possibly different- sets of relevant 

observables with respective expectation values {g a (t)} } {g a (t)}- Let the initial 
state be canonical with respect to the first set {G a }: 

P(t ) = Pc,n[ga(t )] ■ (A.29) 

By assumption, the relevant entropy and the von Neumann entropy are ini- 
tially equal: S re \[g a (t )] = S[p(t )]. Then by unitarity, .^[^(to)] = S'f/?^)], and 
by the general inequality for coarse-graining, .^[^(t)] < S re \[g a (t)]. Hence for 
any later time t > t , and for any choice of the second set of relevant 

observables, it is always 

Sr e i[g a {t)]>S rel [g a {t )] . (A.30) 

This is the second law in its most general form [17]. Whenever both the initial 
and the final macrostate of an isolated system are described by a generalized 
canonical distribution, possibly with respect to different levels of description, 
then the associated final relevant entropy cannot be smaller than the initial 
relevant entropy. In case the evolution of the system leads from one equilibrium 
state to another, the sets {G a } and represent the "old" and "new" 

constants of the motion, respectively. While the second law ensures that the 
final relevant entropy is never smaller than the initial relevant entropy, it does 
not exclude local oscillations of the relevant entropy at intermediate times. 
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Even in the presence of such local oscillations, however, the relevant entropy 
never falls below its initial value. 



The second law follows from reproducibility. As discussed lucidly by 
Jaynes [107], the second law is a consequence of the fundamental requirement 
that experiments on macroscopic systems be reproducible. In such experi- 
ments the microscopic details of the system are generally beyond the control 
of the experimenter; all he can control and observe is some set of macroscopic 
parameters. Hence the preparation of the system amounts to adjusting those 
few expectation values {g a (to)} which characterize the initial macrostate. The 
system is then allowed to evolve -usually until it reaches a new equilibrium- 
and at the end it is again only a limited set of expectation values {g a (t)} } 
namely those which characterize the final macrostate, that will be measured. 
The experiment is reproducible if its outcome (i. e., the final expectation val- 
ues {g a (t)}) is uniquely determined by the initial preparation (i. e., by the 
initial expectation values {(/ a (^o)}); in other words, if the final {g a (t)} can be 
predicted on the basis of the initial {g a (to)}- But a prediction cannot possibly 
contain more information than the data on which it is based. Hence the final 
expectation values cannot carry more information (as to the microstate of the 
system) than do the initial expectation values; which implies that the cor- 
responding final relevant entropy cannot be smaller than the initial relevant 
entropy, Q.E.D. 

A.3.2 H-theorem 

On a Markovian level of description the relevant entropy increases 
monotonically. If and only if the level of description spanned by the rele- 
vant observables {G a } is Markovian, the general argument for the second law 
can be applied iteratively: the expectation values {g a (t + dt)} at time t + dt 
cannot carry more information than the expectation values {g a (t)} at time 
t, based on which they could be predicted. (A closely related argument for 
the _ff-theorem is sketched in figure A.l.) As a result, the associated relevant 
entropy can only increase or stay constant, but never decrease; in short: 



There is a common misconception that the _ff-theorem constitutes a proof 
of the second law. This is wrong. The _ff-theorem is in fact weaker than the 
second law in that it refers to only one set of relevant observables, not two; 
and that it holds only if the level of description is Markovian. On the other 
hand, it is stronger in that it describes the behavior of the relevant entropy 
at all times t, rather than just comparing initial and final values. 



Sr e i[g a (t)] > Vt 
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Pmttl) 



set 



(non-unitary) Markovian 
time evolution 



P(ti)=P m (t 1 ) 



unitary evolution 



A discard 
irrelevant 
information 



P(t 2 ) 



Fig. A.l. One argument for the if-theorem. p m denotes the reduced state associated 
with the Markovian level of description; p denotes the full state. In the Markovian 
limit the diagram must be commutative. Since the lower path involves discarding 
information, entropy must necessarily increase: S m {t2) > S m (ti), Q.E.D. 

Information flows from slow to fast degrees of freedom. The H- 

theorem reveals once more the special status of the Markovian level of de- 
scription: not only is the Markovian level particularly useful for explicit cal- 
culations, it is also distinguished conceptually. As we discussed before, this 
distinction is linked to the existence of widely disparate time scales. On the 
Markovian level of description "relevant" means "slow," while "irrelevant" 
means "fast." The growth of the Markovian relevant entropy indicates that 
the amount of information carried by the slowly varying observables steadily 
decreases. At the same time the constancy of the von Neumann entropy shows 
that complete information about the system is retained in a full microscopic 
description. An obvious interpretation is that in the course of the system's evo- 
lution, information is continuously being transferred from slow to fast degrees 
of freedom. It is this flow of information which is perceived as irreversible. The 
information is not really lost; it only becomes inaccessible to a certain coarse- 
grained (namely the Markovian) level of description. In the example of a dilute 
classical gas, information is being transferred from single-particle observables 
to many-particle correlations. Since in general these correlations will not be 
measured, part of the information about the system becomes experimentally 
inaccessible. 



As the Markovian entropy approaches the equilibrium entropy, the 
system thermalizes. For brevity, let us denote by p m (t) the relevant part of 
the statistical operator which is associated with a Markovian level of descrip- 
tion, and by S m (t) the associated Markovian entropy. The Markovian level 
of description, being spanned by all slowly varying observables, contains as a 
subspace the equilibrium level of description which is spanned by all constants 
of the motion. As a consequence, the Markovian entropy can never be greater 
than the equilibrium entropy: 

S m (t)<S eq W . (A.32) 
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While the Markovian entropy keeps increasing due to the _ff-theorem, the 
equilibrium entropy is a constant. Whether or not the Markovian entropy 
eventually approaches the equilibrium entropy as t — > oo, depends on the 
particular physical process. If it does, i. e., if 

S m (t) -> S eq (A.33) 
and hence 

S{p e(i \p m {t)) -> (A.34) 

then, because the relative entropy is an appropriate distance measure, one 
may conclude that the Markovian distribution approaches the equilibrium 
distribution. The level of description may then be contracted from the original 
Markovian to the equilibrium level. This means that the system's macrostate, 
originally described on the larger Markovian level of description, is now char- 
acterized completely by the equilibrium distribution and hence solely by the 
expectation values of the constants of the motion. This level contraction is 
the mathematical manifestation of the physical "equilibration," or "thermal- 
ization," of the system. 

The above line of argument reveals the importance of the entropy concept even 
away from equilibrium. Being an appropriate distance measure, the entropy 
difference [S eq — S m (t)] is the proper diagnostic tool to describe the approach 
to equilibrium. 

Thermalization may occur in several stages. Between the original (M.^) 
and the equilibrium (£) level of description there may lie several intermediate 
Markovian levels of description, 

M {0) D ...D M (t) D M (t+1) D ...D M (n) D £ , (A.35) 

each spanned by observables which evolve on successively longer time scales. 
Associated with these intermediate levels are reduced states p^(t) and en- 
tropies Sf$(t). The thermalization of the system then occurs in various stages, 
each stage amounting to a level contraction J\4^ — > Again, the en- 

tropy difference 

" S$(t)] = S( P g +1 \t)\pg\t)) (A.36) 

is the proper diagnostic tool to test the feasibility of such a contraction. The 
thermalization is complete as soon as one may contract J\4^ — > £. 
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The various intermediate level contractions occur on successively longer time 
scales; they represent different physical "regimes" of the macroscopic evolu- 
tion. Two examples are the kinetic and hydrodynamic regimes in the evolution 
of a dense plasma. The former is described on the Boltzmann level of descrip- 
tion, spanned by all single-particle observables; whereas the latter is described 
on the hydrodynamic level of description, spanned only by the observables of 
local particle, energy, and momentum density. As the plasma evolves, it first 
reaches local equilibrium, permitting to contract the Boltzmann to the hy- 
drodynamic level of description; then, on a generally much longer time scale, 
it approaches global equilibrium, eventually permitting to contract the hy- 
drodynamic to the equilibrium level of description. These regimes and their 
associated time scales can be identified by studying solely the evolution of the 
relevant entropies Sf$(t); a fact which underlines once more the role of the 
entropy powerful diagnostic tool. 

The i7-theorem is experimentally relevant. One may wonder about 
the experimental relevance of the _ff-theorem, given that the Markovian level 
of description usually differs from the experimental level. However, for an ob- 
servable to be measurable in practice it is usually necessary that it vary slowly. 
Hence the Markovian level typically contains the experimental level (Fig. A. 2). 
For this typical behavior of the associated entropies is shown in Fig. 

A. 3. Both the von Neumann and the equilibrium entropy remain constant 
in time. In contrast, the Markovian entropy increases monotonically due to 
the _ff-theorem. As t — > oo the Markovian entropy need not always approach 
the equilibrium entropy; however, if it does, then the system has thermal- 
ized. About the behavior of the experimental entropy no general statement 
can be made — it may or may not have local extrema. But as long as the 
experimental level of description is contained in the Markovian level, it is 
S m < S exp < S eq . Then thermalization on the Markovian level immediately 
implies thermalization on the experimental level. 

That the experimental entropy may actually decrease temporarily, can be il- 
lustrated with a simple example [36]. Imagine, for instance, two clouds of a 
very dilute, very weakly interacting gas, each initially in local equilibrium, 
moving towards one another with large average velocities. Let us assume that 
only the local particle, energy, and momentum densities are measured experi- 
mentally. On this level of description the associated experimental entropy will 
sharply increase as soon as the two clouds begin to overlap; but after the en- 
counter, it will fall back to practically its initial value. This decrease reflects 
the non-Markovian nature of the experimental level of description: the values 
of the local densities at time t + dt depend not only on their values at time 
t, but also on their past history before the encounter. In contrast, the Boltz- 
mann entropy -which is associated with a Markovian level of description- will 
increase monotonically throughout the evolution of the system. 
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general case 




typical case 



Fig. A. 2. Levels of description of a physical system. 1 denotes the unit operator, 
C the constants of the motion, E the experimentally measured observables, S the 
slow observables, and A all observables. The arrows indicate the flow of information 
in the Markovian limit (if-theorem) . 
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Fig. A. 3. An example for the behavior of various entropies. Both the von Neumann 
and the equilibrium entropy remain constant in time. The Markovian entropy in- 
creases monotonically (if-theorem) . No such general statement can be made about 
the behavior of the experimental entropy; it may or may not have local extrema. 
However, the experimental level of description is typically contained in the Marko- 
vian level, and hence S'exp > S m . 



72 



References 



[1] S. Nakajima. On quantum theory of transport phenomena. Prog. Theor. 
Phys., 20:948, 1958. 

[2] R. Zwanzig. Ensemble method in the theory of irreversibility. J. Chem. Phys., 
33:1338, 1960. 

[3] R. Zwanzig. Memory effects in irreversible thermodynamics. Phys. Rev., 
124:983, 1961. 

[4] R. Zwanzig. Statistical mechanics of irreversibility. In W. E. Brittin, B. W. 
Downs, and J. Downs, editors, Boulder Lectures in Theoretical Physics, Vol. 
3. Interscience, 1961. 

[5] R. Zwanzig. On the identity of three generalized Master equations. Physica 
(Utrecht), 30:1109, 1964. 

[6] R. Zwanzig. Nonlinear generalized Langevin equations. J. Stat. Phys., 9:215, 
1973. 

[7] H. Mori. Transport, collective motion, and Brownian motion. Prog. Theor. 
Phys., 33:423, 1965. 

[8] H. Mori. A continued-fraction representation of the time-correlation functions. 
Prog. Theor. Phys., 34:399, 1965. 

[9] B. Robertson. Equations of motion in nonequilibrium statistical mechanics. 
Phys. Rev., 144:151, 1966. 

[10] B. Robertson. Equations of motion of nuclear magnetism. Phys. Rev., 153:391, 
1967. 

[11] B. Robertson. Equations of motion in nonequilibrium statistical mechanics. 
II. energy transport. Phys. Rev., 160:175, 1967. (Erratum: 166:206, 1968). 

[12] B. Robertson. Quantum statistical mechanical derivation of generalized 
hydrodynamic equations. J. Math. Phys., 11:2482, 1970. 

[13] B. Robertson. In R. D. Levine and M. Tribus, editors, The Maximum Entropy 
Formalism. M.I.T. Press, 1978. 

[14] K.-H. Li. Physics of open systems. Phys. Rep., 134:1, 1986. 

[15] F. Haake. Statistical Treatment of Open Systems by Generalized Master 
Equations. Springer, Berlin, 1973. 

[16] H. Grabert. Projection Operator Techniques in Nonequilibrium Statistical 
Mechanics. Springer, Berlin, Heidelberg, 1982. 

[17] E. Fick and G. Sauermann. The Quantum Statistics of Dynamic Processes. 
Springer, Berlin, Heidelberg, 1990. 



73 



[18] D. N. Zubarev. Boundary conditions for statistical operators in the theory of 
nonequilibrium processes and quasiaverages. Theor. Math. Phys., 3:505, 1970. 

[19] D. N. Zubarev. Nonequilibrium statistical thermodynamics. Plenum, New 
York, 1974. 

[20] G. Ropke. Statistische Mechanik fur das Nichtgleichgewicht. Physik-Verlag, 
Weinheim, 1987. 

[21] D. N. Zubarev and M. V. Tokarchuk. Nonequilibrium thermofield dynamics 
and the nonequilibrium statistical operator method. Theor. Math. Phys., 
88:876, 1991. 

[22] H. Chu and H. Umezawa. A unified formalism of thermal quantum field theory. 
Int. J. Mod. Phys. A, 9:2363, 1994. 

[23] P. A. Henning. Thermo field dynamics for quantum fields with continuous 
mass spectrum. Phys. Rep., 253:235, 1995. 

[24] S. R. de Groot, W. A. van Leeuwen, and C. G. van Weert. Relativistic Kinetic 
Theory. North-Holland, Amsterdam, 1980. 

[25] P. Carruthers and F. Zachariasen. Quantum collision theory with phase-space 
distributions. Rev. Mod. Phys., 55:245, 1983. 

[26] J. Schwinger. Brownian motion of a quantum oscillator. J. Math. Phys., 2:407, 
1961. 

[27] L. P. Kadanoff and G. Baym. Quantum Statistical Mechanics. Benjamin, New 
York, 1962. 

[28] L. V. Keldysh. Diagram technique for nonequilibrium processes. Zh. Eksp. 
Tear. Fiz., 47:1515, 1964. [JETP 20 (1965) 1018]. 

[29] K.-C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu. Equilibrium and nonequilibrium 
formalisms made unified. Phys. Rep., 118:1, 1985. 

[30] N. N. Bogoliubov. In J. de Boer and G. E. Uhlenbeck, editors, Studies in 
Statistical Mechanics, Vol. I. North-Holland, 1962. 

[31] H. Spohn. Kinetic equations from Hamiltonian dynamics: Markovian limits. 
Rev. Mod. Phys., 53:569, 1980. 

[32] I. Prigogine. Non-Equilibrium Statistical Mechanics. Interscience, New York, 
London, 1962. 

[33] I. Prigogine. Introduction to Thermodynamics of Irreversible Processes. 
Interscience, New York, London, 1969. 

[34] E. T. Jaynes. Papers on Probability, Statistics and Statistical Physics. Kluwer 
Academic, Dordrecht, 1989. ed. by R. D. Rosenkrantz. 

[35] W. Brenig. Statistical Theory of Heat. Springer, Berlin, Heidelberg, New York, 
1989. 



74 



[36] R. Balian. From Microphysics to Macrophysics, Vol. II. Springer, Berlin, 
Heidelberg, 1992. 

[37] G. F. Bertsch, H. Kruse, and S. D. Das Gupta. Boltzmann equation for heavy 
ion collisions. Phys. Rev. C, 29:673, 1984. 

[38] P. Danielewicz. Quantum theory of nonequilibrium processes, I. Ann. Phys. 
(N.Y.), 152:239, 1984. 

[39] C. Greiner, K. Wagner, and P.-G. Reinhard. Memory effects in relativistic 
heavy ion collisions. Phys. Rev. C, 49:1693, 1994. 

[40] C. Greiner, K. Wagner, and P.-G. Reinhard. Finite memory in the collision 
processes of afermionic system and its effect on relativistic heavy ion collisions. 
In W. Greiner, H. Stocker, and A. Gallmann, editors, Hot and dense nuclear 
matter. Plenum, 1994. 

[41] J. M. Eisenberg and G. Kalbermann. Pair production in transport equations. 
Phys. Rev. D, 37:1197, 1988. 

[42] I. Bialynicki-Birula, P. Gornicki, and J. Rafelski. Phase-space structure of the 
Dirac vacuum. Phys. Rev. D, 44:1825, 1991. 

[43] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, and E. Mottola. Pair 
production in a strong electric field. Phys. Rev. Lett., 67:2427, 1991. 

[44] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, and E. Mottola. Fermion 
pair production in a strong electric field. Phys. Rev. D, 45:4659, 1992. 

[45] J. Rau. Pair production in the quantum Boltzmann equation. Phys. Rev. D, 
50:6911, 1994. 

[46] B. L. Hu and H. E. Kandrup. Entropy generation in cosmological particle 
creation and interactions: A statistical subdynamics analysis. Phys. Rev. D, 
35:1776, 1987. 

[47] J. Rau. Entropy production via particle production. In J. Letessier, 
H. H. Gutbrod, and J. Rafelski, editors, Hot Hadronic Matter: Theory and 
Experiment. Plenum, 1995. 

[48] J. L. Lebowitz. Boltzmann's entropy and time's arrow. Phys. Today, Sept. 
1993, p. 32ff. See also letters in Phys. Today, Nov. 1994, p. llff. 

[49] H. D. Zeh. Toward a quantum theory of observation. Found. Phys., 3:109, 
1973. 

[50] W. H. Zurek. Pointer basis of quantum apparatus: Into what mixture does 
the wave packet collapse? Phys. Rev. D, 24:1516, 1981. 

[51] W. H. Zurek. Environment-induced superselection rules. Phys. Rev. D, 
26:1862, 1982. 

[52] J. J. Halliwell. Decoherence in quantum cosmology. Phys. Rev. D, 39:2912, 
1989. 



75 



[53] A. Albrecht. Investigating decoherence in a simple system. Phys. Rev. D, 
46:5504, 1992. 

[54] B. L. Hu, J. P. Paz, and Y. Zhang. Quantum Brownian motion in a general 
environment: Exact master equation with nonlocal dissipation and colored 
noise. Phys. Rev. D, 45:2843, 1992. 

[55] J. B. Hartle. Spacetime coarse grainings in nonrelativistic quantum mechanics. 
Phys. Rev. D, 44:3171, 1991. 

[56] M. Gell-Mann and J. B. Hartle. Classical equations for quantum systems. 
Phys. Rev. D, 47:3345, 1993. 

[57] P. Fulde. Electron Correlations in Molecules and Solids. Springer, Berlin, New 
York, 1993. 2nd ed. 

[58] K. Burnett, J. Cooper, R. J. Ballagh, and E. W. Smith. Collisional 
redistribution of radiation. I. The density matrix. Phys. Rev. A, 22:2005, 
1980. 

[59] L. P. Csernai. Introduction to Relativistic Heavy Ion Collisions. Wiley, 
Chichester, 1994. 

[60] R. Kubo, M. Toda, and N. Hashitsume. Statistical Physics II. Springer, Berlin, 
Heidelberg, 1985. 

[61] R. Kubo. Statistical-mechanical theory of irreversible processes. I. J. Phys. 
Soc. Japan, 12:570, 1957. 

[62] R. Kubo. Lectures in Theoretical Physics I. Interscience, New York, 1959. 

[63] K. Kawasaki and J. D. Gunton. Theory of nonlinear transport processes: 
Nonlinear shear viscosity and normal stress effects. Phys. Rev. A, 8:2048, 
1973. 

[64] P. A. M. Dirac. Note on exchange phenomena in the Thomas atom. Proc. 
Camb. Phil. Soc, 26:376, 1930. 

[65] R. Balescu. Equilibrium and nonequilibrium statistical mechanics. Wiley, New 
York, 1975. 

[66] R. D. Richtmyer. Principles of Advanced Mathematical Physics, Vol. I. 
Springer, New York, 1978. 

[67] A. L. Fetter and J. D. Walecka. Quantum Theory of Many-Particle Systems. 
McGraw-Hill, New York, 1971. 

[68] E. Fermi. Nuclear Physics. Univ. Chicago Press, Chicago, 1950. 

[69] L. Van Hove. Quantum-mechanical perturbations giving rise to a statistical 
transport equation. Physica (Utrecht), 21:517, 1955. 

[70] L. Van Hove. The approach to equilibrium in quantum statistics. Physica 
(Utrecht), 23:441, 1957. 



76 



[71] P. Drude. Zur Elektronentheorie der Metalle. Ann. Phys. (Leipzig), 1:566, 
1900. 

[72] P. Drude. Zur Elektronentheorie der Metalle. Ann. Phys. (Leipzig), 3:369, 
1900. 

[73] N. W. Ashcroft and N. D. Mermin. Solid State Physics. Saunders College, 
Philadelphia, 1976. 

[74] F. Sauter. Uber das Verhalten eines Elektrons im homogenen elektrischen 
Feld nach der relativistischen Theorie Diracs. Z. Phys., 69:742, 1931. 

[75] W. Heisenberg and H. Euler. Folgerungen aus der Diracschen Theorie des 
Positrons. Z. Phys., 98:714, 1936. 

[76] J. Schwinger. On gauge invariance and vacuum polarization. Phys. Rev., 
82:664, 1951. 

[77] A. I. Nikishov. Pair production by a constant external field. Sov. Phys. JETP, 
30:660, 1970. 

[78] M. S. Marinov and V. S. Popov. Pair production in an electromagnetic field 
(case of arbitrary spin). Sov. J. Nucl. Phys., 15:702, 1972. 

[79] V. S. Popov. Pair production in a variable and homogeneous field as an 
oscillator problem. Sov. Phys. JETP, 35:659, 1972. 

[80] A. M. Perelomov. Some remarks on fermion pair creation in an alternating 
external field. Phys. Lett. A, 39:353, 1972. 

[81] A. M. Perelomov. Creation of fermion pairs in a homogeneous alternating 
external field. Theor. Math. Phys., 19:368, 1974. 

[82] M. S. Marinov and V. S. Popov. Electron-positron pair creation from vacuum 
induced by variable electric field. Fortsch. Phys., 25:373, 1977. 

[83] J. Rau. The apparent irreversibility of particle creation: A study of time scales 
and of the mechanisms responsible for entropy production in quantum field 
theory. PhD thesis, Duke University, 1993. 

[84] N. N. Bogoliubov. A new method in the theory of superconductivity. III. Zh. 
Eksp. Teor. Fiz., 34:73, 1958. [JETP 7 (1958) 51]. 

[85] N. D. Birrell and P. C. W. Davies. Quantum fields in curved space. Cambridge 
University Press, Cambridge, 1984. 

[86] C. Zener. Non-adiabatic crossing of energy levels. Proc. Roy. Soc. A, 137:696, 
1932. 

[87] A. M. Dykhne. Adiabatic perturbation of discrete spectrum states. Sov. Phys. 
JETP, 14:941, 1962. 

[88] L. D. Landau and E. M. Lifshitz. Quantum Mechanics. Pergamon, Oxford, 
1965. 2nd ed. 



77 



[89] G. A. Hagedorn. Proof of the Landau-Zener formula in an adiabatic limit with 
small eigenvalue gaps. Commun. Math. Phys., 136:433, 1991. 

[90] A. Casher, H. Neuberger, and S. Nussinov. Chromoelectric-flux-tube model 
of particle production. Phys. Rev. D, 20:179, 1979. 

[91] V. S. Olkhovsky and E. Ricami. Recent developments in the time analysis of 
tunnelling processes. Phys. Rep., 214:339, 1992. (The proper definition of a 
"tunnelling time" is still controversial.). 

[92] F. Cooper, J. M. Eisenberg, Y. Kluger, E. Mottola, and B. Svetitsky. Particle 
production in the central rapidity region. Phys. Rev. D, 48:190, 1993. 

[93] N. N. Bogoliubov. Kinetic equations. J. Phys. (USSR), 10:265, 1946. 

[94] C. E. Shannon. Bell System Tech. J., 27:379,623, 1948. 

[95] C. E. Shannon and W. Weaver. The Mathematical Theory of Communication. 
University of Illinois Press, Urbana, 1949. 

[96] A. Renyi. Wahrscheinlichkeitsrechnung. Deutscher Verlag Wiss., Berlin, 1979. 
6. Aufl. 

[97] S. Kullback. Information Theory and Statistics. John Wiley, New York, 1959. 

[98] E. T. Jaynes. Information theory and statistical mechanics. In K. W. Ford, 
editor, Statistical Physics. Benjamin, 1963. (1962 Brandeis Theoretical Physics 
Lectures). 

[99] W. Thirring. A Course in Mathematical Physics 4: Quantum Mechanics of 
Large Systems. Springer, New York, Wien, 1983. 

[100] E. T. Jaynes. Where do we stand on maximum entropy? In R. D. Levine and 
M. Tribus, editors, The Maximum Entropy Formalism. M.I.T. Press, 1978. 

[101] E. T. Jaynes. Information theory and statistical mechanics. Phys. Rev., 
106:620, 1957. 

[102] E. T. Jaynes. Information theory and statistical mechanics. II. Phys. Rev., 
108:171, 1957. 

[103] A. Katz. Principles of Statistical Mechanics. Freeman, San Francisco, London, 
1967. 

[104] R. Balian and N. L. Balazs. Equiprobability, inference, and entropy in 
quantum theory. Ann. Phys., 179:97, 1987. 

[105] E. T. Jaynes. Prior probabilities. IEEE Trans. Systems Sci. Cyb., 4:227, 1968. 

[106] E. T. Jaynes. The well-posed problem. Found. Phys., 3:477, 1973. 

[107] E. T. Jaynes. Gibbs vs Boltzmann entropies. Am. J. Phys., 33:391, 1965. 



78 



